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ABSTRACT 



We report on the A = 5 — 36/iin Spitzer Infrared Spectrograph spectra of 79 young 
stellar objects in the very young nearby cluster NGC 1333. NGC 1333's youth enables 
the study of early protoplanetary disk properties, such as the degree of settling as 
well as the formation of gaps and clearings. We construct spectral energy distributions 
(SEDs) using our IRS data as well as published photometry and classify our sample into 
SED classes. Using "extinction-free" spectral indices, we determine whether the disk, 
M-f envelope, or photosphere dominates the spectrum. We analyze the dereddened spectra 

of objects which show disk dominated emission using spectral indices and properties 
of silicate features in order to study the vertical and radial structure of protoplanetary 
disks in NGC 1333. At least nine objects in our sample of NGC 1333 show signs of large 
(several AU) radial gaps or clearings in their inner disk. Disks with radial gaps in NGC 
1333 show more-nearly pristine silicate dust than their radially continuous counterparts. 
We compare properties of disks in NGC 1333 to those in three other well studied regions, 
Taurus-Auriga, Ophiuchus and Chamaeleon I, and find no difference in their degree of 
CH \ sedimentation and dust processing. 

Subject headings: accretion disks - circumstellar matter - dust, extinction - infrared: 
stars - open clusters and associations: individual (NGC 1333, Taurus-Auriga, Ophi- 
uchus, Chamaeleon I) - planetary systems: formation, protoplanetary disks — stars: 
formation 
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Introduction 



NGC 1333 was first identified as an optical reflection nebula at the western end of the Perseus 
molecular cloud. The moniker NGC 1333 is also used to identify the molecular cloud an d young stel- 



lar cluster associated with the optical reflection nebula along w ith the dark cloud L1450 (JHerbig fc Kameswara Rao 
1972J : lrIerbiglll974! : lstrom et al.lll97d : lwalawender et alihoosh . NGC 1 333 is one of the best studied 



very 



1996 



young; clusters of acti ve low- to intermediate- mass star formation (lAspin et al.lll994l : lLada et al 



Wilking et al.l 12004 ). Imaging surveys in the near- and mid-infrared have r e vealed NGC 1333 



contains both a norther n and southern young stellar cluster (jAspin et al.l LL994J; lLada et al.l Il99d ; 
Gutermuth et al.ll2008bl ). Previous studies have also shown NGC 1333 is heavily extinguished and 
contains a substanti al fraction of deep l y embedded young ste l lar objects (YSOs) a s well as many 
prominent outflo ws (lLada et al.l Il99a : iKnee fc Sandelll 120001 ; IWinston et al.1 120091). The Spit zer 



Space Telescope ( Werner et al.1 12004 ) infrared view of NGC 1333 in iGutermuth et al.l ( 2008bl ) is 



streaked with emission at 4.5/im, which is dominated by H2 lines arising in shocks from YSO 
outflows colliding with ambient cloud material. These energetic outflows originate from actively 
accreting protostars with infalling envelopes (Class 0, I). As they are surrounded by dusty en- 
velopes, these protostars as well as more evolved pre-main-sequence stars with circumstellar disks 
(Class II) emit in the mid-infrared and are well detected by Spitzer. Pre-main-sequence stars with 
little to no dusty disk or envelope (Class III) are n ot as easily detected in the mid-infrared and hav e 



been found in NGC 1333 through X-ray surveys (JGetman et al. 



2002 



Winston et al.ll2009l . l2010h . 



Many previous studies have addres sed the age of NGC 1 333. In a near-infrared (NIR) spectro- 
graphs and X-ray study of NGC 1333, IWinston et al.l (J2009I ) estimated most objects in NGC 1333 



to be <3 Myr, and found an age spread from 1 to 10 Myr b a sed on isochrones by iBaraffe et al. 
(|1998l ). An optical and NIR study of NG C 1333 by lAspinl ()2003l ) found similar results using 



isochrones by lD'Antona fe Mazzitellil f 19971 ) namely an age spread from 7xl0 4 yr to 10 Myr with 
many sources around 3 Myr. They interpreted this spread in ages to be evidence of multiple epochs 
of star formation potentially spurred by turbulence generated by the many prominent outflows in 
NGC 1333. Oth ers have also sugge sted turbulence driven star formation could have taken place in 
NGC 1333 (e.g. buillen et allbood ). Studies of the low mass (0.25 - 0.02 M ) population of NGC 
1333 esti mate a median age of ^ 0- 3 My r for that population using the low mass stellar evolution 



tracks by iD'Antona fc Mazzitellil (J1997J ) (jGreissl et al. 



2007 



Wilking et al.l 12004 ) . However, for 



objects less than 1 Myr old age estimates based on isochrones are more uncertain than for older 
objec ts as the location of the birth line (t = 0) is ambiguous (JD'Antona &: Mazzitellil 119971 : ISiess 
200 ll ) . For t his reason, class fractions are a common tool used to get an idea of the relative youth 
of a cluster. IGutermuth et al.l (|2008bl ) show NGC 1333 is very young due to its large disk fraction 
(83%±11%) which they define as the fraction of total objects which s how a large infrar ed excess 
at certain wavelengths indicative of emission from a circumstellar disk. lLada et al.l (119961 ) estimate 
the age of NGC 1333 to be very young based on the similarity of its i<f-band luminosity function 
to that of the Trapezium (< 1 Myr), as well as its large fraction of infrared sources (61%) which 
suggest an age of < 1-2 Myr. The protostar and Class II populations also have the same peak near- 



est neighbor distance (0.045 pc), m eaning stars do not appe ar to have moved far from where they 
formed, another indication of youth (JGutermuth et al.ll2008bl ). By these measures, NGC 1333 is one 
of the youngest clusters observed comprehensively by Spitzer. Based on disk fraction, NGC 1333 is 
younger than three we l l-stud ied star forming regions a lso observed by Spitze r, namely Taurus (~ 
1-2 Mvr: IFurlan et al.l ( 2006)). Ophiuchus (~ 0.8 Mvr: iMcClure et all (|201oh ). and Chamaeleon I 
(~ 2 Myr; iManoj et al.l ()201ll )). We adopt an age of < 1 Myr for this study. 



The age of a cluster, while difficult to determine, is necessary to understand the timescales of 
star formation and disk evolution processes. Observations show disk dissipation is well underway in 



3-5 Myr, as on this timescale infrared excess around young stellar objects drops by half (Ha isch et al 



2001 



Muzerolle et al. 



Hernandez et al.ll2008l ) and accretion rate falls below 10 10 M Q yr x for low mass stars (e.g. 



2003). Therefore the study of regions significantly younger than this allow for in- 



vestigation of early disk evolution and dissipation processes and can further constrain the timescales 
of these phenomena. Circumstellar dust and gas predominately emit in the infrar ed making it the 
ideal wavelength range to investigate processes in disk evolution. Models predict fWeide nschilling 



1997 



Harker &: Desch 



servations confirm (Furlan et al. 



2009 



200 



2: 



Gai 



Bouwman et al. 



2004 : iDullemond fc Domin 



2005; 



Kessler-Silacci et al. 



20081 : ISargent et al. 



2009 



ik 



200 



5; 



D'Alessio et al. 



2006: ISicilia-Aguilar et al. 



Olofsson et al.li2009l . l2010l ) that grain growth, dust 



)6) 



2006) and ob- 



2007, 



2008 



processing and dust settling in disks proceed on timescales <1 Myr. 

Giant planets are predicted to form within a few million years throug h collisional growth 



of grains into multi-M m cores which may then ac crete gas from the disk (jPollack et al. 



Weidenschilling 



2008 



Dodson-Robinson et al. 



200 



1996 



They may also form th r ough gravit ational 

The 



2001 



instability of gas and dust in the disk in <10 5 yr (JGoldreich fc Ward! 1 19731 : iBosa 
gravitational interaction between a disk and a giant planet would clear out mate rial in its or- 
bit thereby opening a gap or forming an inner clearing (e.g. iBorderies et al.lll989l ). Disks with 



radial gaps and holes have been detected mostly via deficit s in their spectra l energy distribu 



tions (SEDs) compa r ed to the median for their population (IStrom et al. 



D'Alessio et al. 
20091 : 



Merin et al.l 



2005: . 
|2010J 



Brown et al.l 120071 : ICalvet et al.l J2005I : JEspaillat et al 



1989; 



Rice et al. 



2007 



2008 



2003; 



Kim et al 



I) , but in several cases h ave been confirmed by mm- wave interferometric imag- 



ing (e.g. lPietu et al.ll2006l : iHughes et al.ll2009l ). The leading explanation for such gaps and clearings 
in disks is the gravitational and tidal effect of companions with the centr al star. In the paradigm 



"central clearing" disk around CoKu Tau/4 the com panion is a star (jlreland &: 



but in most others it is likely to be a gas-giant planet (JArtymowicz k, Lubow 



Krau 



1994 



200. 



Calvet et al 



2005l ; lKim et al.ll2009l ; iPott et al.ll2010i ). However, several disk-dissipation mechanisms which clear 
disks from the inside out have been invoked to exp lain central clearing s, among them photoe- 



vaporation (IClarke et al.ll200ll : lAlexander et al.l 120061 ) and MRI draining (jBalbus fc Hawlevl Il991 



Chiang fc Murrav-Clavl 120071 ) . These methods generally fail to account for gaps separating in 



ner and outer disks , which outnumber the central clearings (IFurlan et al.l 120091 : iKim et al. 



McClure et al. 



2010 



Manoj et al. 



2011 



2009 



Disks with holes or gaps are historically called transi- 
tional disks (TDs) as they appear to be in transition between a star with a full disk of dust and 



gas arou nd it (Class II object) and a star where l ittle to no emission is s een from a disk(Class III 
object) (jStrom et al.l Il989l ; ISkrutskie et al.l ll99Q ; iGauvin fc Stromlll992r). Disks with radial gaps 
are a type of t r ansitional disk known as 'pre-transitional disks' (PTDs; lEspaillat et al.l 120071 . 120081 ; 
Espaillatl |2009j; lEspaillat et al.l |2010| ) . The interaction of one or more planets with the disk is key 
to understanding the timescales and processes of planet formation as well as disk dissipation. 

Here we present the spectra of 81 objects in NGC 1333. In § [2] we discuss the observations, 
selection, completeness, and data reduction for our sample. We correct our YSO sample for inter- 
stellar extinction in § [3] and using their observed spectral indices we classify our sample into SED 
classes and evolutionary states in § [H In § [5] we characterize the settling in different regions of 
our disks as well as identify objects with interrupted radial disk structure. We also compare the 
distribution of various spectral properties and the median spectrum of disks in NGC 1333 to those 
in Taurus, Ophiuchus and Chamaeleon I. Further discussion of the classification of TDs and PTDs 
in NGC 1333 is found in § [6j Finally, in §[7] we summarize our results and state our conclusions. 



Observations and data reduction 



2.1. IRS observations 



As part of a Spitzer GTO program (PID 40525) we observed 78 objects in NGC 1333 on 2007 
October 5-6 k 9-10, 2008 February 21 to March 2 and 2008 March 30, in IRS campaigns 44, 48 and 
49 respectively. Our sample of YSOs was chosen based on their excess in photometric observations 
from both the Spitzer Infrared Ar ray Camera (IRA C; iFazio et al,ll2004l ) 8.0/j,m band and Multi 



band Im aging Photometer (MIPS: lRieke et al.ll2004l ) 24.0/um band, published in lGutermuth et al. 



(J2008bl ). We observed all objects with detections {a < 0.2 mag) in the IRAC 8.0/im band and 
MIPS 24^m band which are brighter than 11.54 mag at 8.0/um and 8.5 mag at 2 4 um. For this 



reason we adopt their identifiers for our sample (these identifiers were also used in lWinston et al. 



2009 



2010). The basic data for this sample, namely, [GMM08] SED index, RA, Dec, AOR ID and 



alternate names, are listed in Table [TJ 



This selection excluded a total of 55 objects from the iGutermuth et al.l (J2008bl ) YSO sample: 
39 had no MIPS 24//m photo metry, 2 had no IRAC 8. 0iim. photometry and 14 were too faint at 
either 8.0/im or 24/um. Of the IGutermuth et al.l (j2008bl ) sample, 4 objects were already scheduled 
to be observed with IRS and therefore these objects were not observed in program 40525. One of 
these conflicts was the class protostar, IRAS 2A, which was observed in Spitzer program 2. The 
other three, ASR 118, SSTc2d J032924.1+311958, and SSTc2d J032929.3+311835, were observed in 
program 30843 on 2006 September 18, 2007 March 9 and 21. These objects wer e selected from thei r 
photometry as possible 'cold disks' which have inner holes and a r e prese nted in lMerfn et al.l (|2010l ) 
indexed there by #2, #4 and #5 respectively. IGutermuth et al.l (J2008bl ) identifies these objects as 
#67, #122 and #137 an d we include these three observations with IRS in our sample using the 



Gutermuth et al. 



(J2008bl ) identifiers. This addition brings the total number of observations to 81. 



Selecting for objects with excess in the 8.0/um IRAC and 24/um MIPS bands, biases our sample 
toward e arlier star formati o n stage s (Class 0/1 and II) where massive dusty envelopes and disks are 
present. iGutermuth et al.l (|2008bl ) find the 137 cluster members identified with IRAC and MIPS 
to represent 83%±11% of the total cluster membership. Our sample consists of over half of these 
sources (N = 81) which is ~49% of the extrapolated total cluster membership. 

We observed each object from 5 — 38/um using the Staring Mode of the Spitzer IRS in both 
orders of the Short-Low (SL; 5.2 — 14.5/zm) and Long-Low (LL; 14 — 38/^m) resolution modules 
(A/AA ~ 60 — 120). From the Spitzer Science Center (SSC)'s basic calibrated data (BCD) products 
using pipeline S18.7, we extracted a spe ctrum using the IR S instrument team's Spectral Modeling, 
Analysis and Reduction Tool (SMART; iHigdon et al.ll2004l ). We corrected for bad or 'rogue' pixels 
by interpolating the flux from two adjacent good pixels in the dispersion direction. 

For each observation, we took two images with the object at different positions, or nods, in 
the slit, each one third of the width of the slit away from the edge. For most images in our s ample , 
we followed the s t andard procedu r e for e xtracting spectra as is done in iFurlan et al.l (|2006l . 120081 ); 
Kim et al.l (|2009l ); iMcClure et al.l ()2010j). We removed background emission by subtracting the 
image with the source at one nod from the image with the source at the opposite nod. We then 
extracted the spectrum from the sky-subtracted 2D image using a column which tapers as the width 
of the point-spread function does. Spectral extraction and removal of background emission for 60 
of these objects proved troublesome because of a combination of bright cloud emission, multiple 
sources in the slit and/or the extinguished, faint nature of the source. For these observations 
with particularly difficult background emission, we individually inspected and carefully extracted 
a spectrum using the SMART optimal extraction tool, AdOpt (JLebouteiller et alJl201Q ). AdOpt 
allows the user to narrow the extraction window as well as fit background emission in that window 
with an n-th order polynomial. Using AdOpt, we fit the background emission in most cases with a 
polynomial of 0th or 1st order, though a few cases required as high as 3rd order. 

With AdOpt, tapered column extraction or a combination of the two methods, we extracted a 
spectrum for each object from each module, order and nod. To calibrate the flux of the spectra, we 
cre ated relative spect ral response functions (RSRFs) by dividing a template spectrum, a Lac (Al 
V; ICohenet allbooah for SL (5.2-14.5/xm) and £Dra (K2 III) for LL (14.0-36/mi) by its tapered 
column or AdOpt extracted IRS spectrum. We multiplied the target spectra by the similarly 
prepared RSRF for each module, order and nod to calibrate the flux from the IRS. For each nod, 
we combined the calibrated spectra from each module and order. To produce a final spectrum, we 
averaged these spectra from the two nods and calculated uncertainty as half of the nod difference. 
Occasionally, slight mispointing resulted in a decrease in flux in the modules which were mispointed, 
causing a mismatch between the SL and LL spectra. We corrected for this by scaling the fainter 
module, which presumably lacks flux due to mispointing, up by a scalar factor. A typical scaling 
factor was 1.1, with a maximum value of 3.5 (#34). As an additional check on our reduction, 
we extracted a spectrum for each object with the optimal extraction software, UR Optimal (W. 
Forrest and C. Tayrien, personal communication) which extracts a spectrum using data from 5, 



7 or 9 pixels around the nod position. The spectrum extracted with tapered column extraction 
and/or AdOpt had smaller nod-difference uncertainties than that of UR Optimal, therefore we use 
the tapered column and/or AdOpt extracted spectrum for our analysis. We estimate the absolute 
photometric uncertainty in flux in the final spectra to be ~5%. 



2.2. Ancillary Data 

We supplement the IRS spectra with broadband photo metry from the litera ture in the J, H 
and K s bands of the Two Micron All Sky Survey (2MASS: Iskrutskie et aDbood ). NGC 1333 was 
also observed as part of the 2MASS extended mission, 2MASS 6X. Data from this mission were 
taken from observations with six times the integration time of standard 2MASS observations. This 
resulted in ~ 1 mag improvement in sensitivity li mit and reduced the u ncertainty for faint objects 
near the detection limit of regular 2MASS data ( Carpenter et al.ll200ll ). Where available, we use 
the 2MASS 6X data over the standard 2MASS from the lit erature. We a lso su pplement our data, 
where available with R and I broadband photometry from lAspin et al.l (|1994l ). We used 2MASS 
K s broadband photometry in our S EDs and whe r e 2MA SS K„-haxv& photome try was not available 
we used K-hsnd photometry from lAspin et al.1 (119941 ) or lLada et al.1 (119961 ). Every object was 
obser yed photometrica ll y with Spitzer and we supplement our data with broadband photometry 
from butermuth et al.l (ho08bl ) in the IRAC 3.6, 4.5, 5.8 and 8.0/xm bands as well as the MIPS 
24//m band. 



2.2.1. Notes on Individual Objects 

Spatially variable, bright sky emission resulted in careful reduction of 60 of the 81 images using 
off nod and off order sky subtraction as well as polynomial fit sky subtraction with AdOpt. Below 
are the details of particularly troublesome objects. We omit objects from our sample which have a 
spectrum that is not representative of a YSO, namely #8 and #38. 

#5: Bright emission lines in the SL spectrum could not be removed through various sky 
subtraction methods; this results in a seemingly discontinuous spectrum. 



#8: Identified as IRAS 4B in iGutermuth et al.1 (J2008bl ). #8 is not identical with the Class 
object but corresponds instead to one of IRAS 4B's outflow lobes. The observation for #8 was 
pointed at the bright emission originating from one of IRAS 4B's outflow cavities to the south 
of the embedded object's sub-mm coordinates. The SED is shown in Figure Q] and it is left out 
of our sample. The SL spectrum is dominated by the emission from the outflow as evidenced by 
the collisionally excited molecular hydrogen emission lines. Th e LL spectrum con tains a forest of 
attenuated, blended water emission lines from IRAS 4B itself (jWatson et al.l 120071 ) though it lies 
about 5" away, at the edge of the slit. 



#23: Careful extraction was unable to remove H2 and polycyclic aromatic hydrocarbon (PAH) 
emission features which do not originate from #23. 

#34'- The SL spectrum was very difficult to disentangle from the sky emission and had to be 
scaled up by 3.5 to match the LL spectrum. Also, the SL spectrum possibly suffers from an offset 
in the dispersion direction. 

#37: The LL spectrum is dominated by a bright source near the target coordinates. The flux 
contribution of this nearby bright source dominated the spectrum and could not be removed. For 
this reason we only present the SL spectrum of #37 and only use the SL spectrum in our analysis. 

#38: Based on the presence of redshifted PAH emission #38 is a galaxy. Using the redshifted 
6.2, 7.7, 8.6 and 11.2/um PAH emission features we calculate a redshift from each feature and obtain 
an average redshift of 0.29. The SED of the galaxy labeled here as #38 is shown in Figure [2] and 
it is left out of our sample. 

#52: The 2nd module of SL nod 1 had data missing shortward of 5.43/um as a result of 
the SSC pipeline trying to correct for rows in the detector which were altered by peak-up camera 
saturation. This resulted in an unreliable nod 1 spectrum and for this reason we use only the nod 
2 spectrum for SL2. 

#58: Two close sources in the LL slit provided some difficultly. Narrow window extraction 
at nod position with AdOpt yielded the best result. LL was scaled by 0.7 to match SL. We include 
#58 in our analysis with the caveat that the LL spectrum is likely not solely representative of #58. 

#63: Both nods of the 1st module of SL displayed a striping effect in the dispersion direc- 
tion, deemed 'jail-bars,' likely due to correction for a bright source in the peak-up array. Careful 
extraction using AdOpt was able to recover the source emission and we use this spectrum in our 
analysis. This effect as well as the faint nature of the object is the likely cause of the noise in the 
spectrum of #63. 

#67: Both nods of the 1st module of LL were severely affected by stray emission in the 
detector from a very bright source just outside the LL slit; f or this reason we e xclude the LL1 



data from the SED. This effect is exemplified in its spectrum in lMerm et al.l (120101') as a wide fals e 



30//m feature, also seen in the IRS spectrum of Hn 10 E in Chamaeleon I (JManoi et al.l 12011 



Merin et al.1 (120101 ) classify this source (called therein ASR118) as a 'cold disk,' a term equivalent 
to our TD. However, due to lack of reliable data in LL1, we exclude #67 from our analysis and 
cannot definitively address its radial structure. 

#68: The emission lines in the SL module at 5.5, 6.1, 6.9, 8.0 and 9.6/im are highly spatially- 
variable molecular hydrogen lines from the sky, which could not be removed in the reduction 
process. 

#75, #94, #110 & #133: The spectrum of each of these objects contains one or more spikes 
due to difficult to remove background emission and/or possible remnant bad pixels. 



#84- Sky subtraction in LL proved difficult, the strong dip around 30/iin is not a real feature 
but could not be removed. 

#92 & #102: Variable, bright sky emission and the faint nature of the source cause these 
carefully extracted spectra to be noisy and uncertain. 



3. Extinction Correction 



As discussed in §Q3 NGC 1333 is heavily reddened; therefore in order to understand the intrinsic 
properties of this sample we determined the amount of reddening toward each object in a uniform 
way. To determine reddening, or extinction, in terms of Ay, we utilized 2MASS photometry because 
it is the most complete set of NIR photometry available for this sample. All but 12 objects have 
good detections (quality flags A, B or C) in J, H and K s . We a dopt a method to determine the 
extinction toward these objects similar to that by iMcClurd (|2009l ) . First, we plotted (Figure [3|) the 
observed J — H and H — K s colors of o ur sample against each other and co mpared their distribution 
to that of t he classical T Tauri locus (JMever. Calvet &; Hillenbrandlll997l ) and main sequence stars 
and dwarfs (IKenyon Sz Hartmannlll995l ). Most of our sample's colors are consistent with a reddened 
classical T Tauri star. However, there are a few outliers: 10 objects are on the right and 3 objects 
fall to the left of the reddened T Tauri locus. The colors of the objects to the right of the classical 
T Tauri locus could also be more heavily reddened due to local extinction as would be the case for 
Class 0/1 objects which show envelope dominated emission. 

To determine the extinction towards these objects we first compiled a list of spectral types 
from the literature. The large extinction towards NGC 1333 (Ay > 10) makes optical and near 
infrared spectra and subsequent spectral type difficult to obtain. From a literature search, 43 
objects have known spectral types which are listed in Table [2 Figure 2] shows the distribution of 
known spectral types. For objects with known spectral t ypes we adopted the intrins ic J — H and 
H — K photospheric colors for a given spectral type from IKenyon &: Hartmannl (|1995l ) transformed 
to the 2MASS system to determine the color excess. As we assume photospheric colors for all 
objects with known spectral type, extinctions calculated will be overestimates for objects with 
H or K s band excess. For objects with no spectral type information, we made the assumption 
that all the extinction is foreground extinction and determined the color excesses, E(H — K s ) and 

rpf tj is \ 

E(J — H), by extrapolating the observed color of each object along the reddening vector e(j-H) 
= 1.63, until their intrinsic J — H color and H — K s color fell along the classical T Tauri locus 



(JMever. Calvet &: Hillenbrandll 19971 ). Using the J — H and H — K s color excess, generally expressed 
as E(\\ — A2), we calculated two values of Ay for each object using 



A 



Ay 

A Xo 



V 



x [E(X 1 - A 2 )] 



(1) 



with a value of 



— from the Ry = 5 version of the Mathis ( 1990 ) extinction law with -r^— 



7.75. 



The Ry = 5 extinction l aw is appropri ate for dense molecular clouds (|Mathislll99Cll ). This value of 
Ry is in agreement with ICernisI (|199dl )'s findings for the L1450 dark cloud in NGC 1333. Those 
authors estimate to have ifV~5 toward the center, around SVS 3, and Ry^A elsewhere. 

We adopt either the J — H or H — K s derived Ay for our analysis. For objects with known 
spectral type, the J — H and H — K s derived Ay were generally within 2 mag with a maximum 
difference in one case of 7 mag. For objects with no spectral type, typical differences between the 
two vales of Ay were around 1 mag, with a maximum difference of 9 mag. When the Ay values 
were within 1 mag, we chose the smallest of the two for our analysis. For differences in derived 
Ay greater than 2 mag, we compared the dereddened photometry to the appropriate photosphere 
scaled to J band, for a given spectral type where available, and chose the Ay which gave the most 
reasonable physical excess. In most cases the Ay derived from the J — H band excess was used. 
For the 8 objects without spectral type information and difference in derived Ay greater than 2 
mag, we compared the spectrum and photometry dereddened with the each Ay value and chose 
the value which minimized the presence of silicate and ice absorption features. Two objects had 
a negative Ay based on J — H or H — K s and two objects had a negative Ay based on J — H 
and H — K s ; for these objects we adopted an Ay of 0. For the objects to the left and right of the 
classical T Tauri locus, the J — H and H — K s derived Ay tended to differ by 15 mag and in one 
case by 40 mag; it should be noted there is greater u ncert ainty in their Ay. The Ay value for each 
object, derived from the Ry = 5 version of lMathid (jl990l ) extinction curve, as well as the method 
with which they were derived is noted in columns 2 and 3 of Table [2] and a histogram showing the 
distribution of Ay for our sample is plotted in Figure [5j The objects with no 2MASS J, H or K s 
band measurements are noted with a null value. The median Ay for NGC 1333 is 5.8 mag. 

Overall, thes e values of Ay differ on average about 2 mag from the equivalent values of Ak in 
Gutermuth et al.l ()2008bl ). In only one case do the two values differ by as much as 12 mag. This 
object, #25, falls to the right of the classical T T a uri loc us and has large extinction. In order to 
characterize low level extinction, I Gutermuth et al.l ()2008bl ) derived an Ak from an extinction map 
generated for the 2MAS S H — K s colors of field stars in the region using the method described in 
Gutermuth et al.1 (|2005l ). Because of the high galactic latitude of NGC 1333, extinction maps based 
on field stars are low resolution, in this case resolution ranges from 70" to 200". We attribute the 
discrepancies in our extinction values to this difference in methodology. 

Of the 66 objects with reliable J — H and H — K s broadband colors, 4 9 hav e an Ay > 3, where 



the molecular cloud extincti on law is known to diverge from the M athis (|1990l ) Ry=5 law for the 



diffuse interstellar medium CMcClure 



2009 



Chapman et al J 120091). F or this reason, we deredden 



the broadband photometr y and the IRS spectra with iMcClurd (120091 ) 's composite extinction law, 
which follows the Mathis (|1990l ) Ry=5 law for Ay < 3 and corrects objects with Ay > 3 with 
two different curves, one for 3 < Ay < 8 and one for Ay > 8. This composite extinction law has 
also been used to de redden the IRS spe ctra of young stars in the nearby clusters of Ophiuchus 



(JMcClure et al.ll20ld ) and Chamaeleon I (JManoj et all 12011 
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4. Characterization of the Sample 

4.1. SED Classification 

Classification of a YSO is typ ically accomplished using the slope of its SED spanning from the 
near- to mid- infrared (|Ladalll987l ). We find the logarithmic slope, n, or spectral index, in a given 
wavelength range, henceforth expressed as n\ 1 -\ 2 , which is given as 



n \i-\ 2 



log(A 2 F A2 )-log(AiF Ai ; 
log(A 2 ) - log(Ai) 



(2) 



An IR slope, typically from 2 to 25/jm, greater than 0.3 corresponds to that of a Class /object, 
"-2-25 between -0.3 and 0.3 is that of Flat Spectrum object, 77,2-25 between -1.6 and -0.3 defines a 
Class //object and 71-2-25 l ess than -1.6 classifies a Class ///object. The empirical classification 



of th e SED of YSOs has become inextricably tied with an evolutionary sequence (j Adams et al. 
19871 ) as discussed in § [TJ therefore proper classification is important in order to understand the 
evolutionary state of these systems. 

We empirically classify the observed SEDs of the objects in our sample with 2MASS K s band 
photometry using 77-2-25- The flux at 2/xm is the 2MASS K s photometric flux (A=2.17/im) and 
the flux at 25/iin is the integrated flux of the observed IRS spectrum from 23.5-26.5/i-m. SED 
classification from 2 to 25/jm could not be done for 8objects because they lack detection in the K s 
band of 2MASS and for 1 object (#67) because of unreliable data near 25/_tm. We find 16 Class 
I objects, 14 Flat Spectrum (henceforth, FS) objects and 40 Class II objects based on 712-25 SED 
classification. No objects were classified as Class III based on their 712-25- This is not surprising 
as Class III objects would not be bright in the mid-infrared and we selected our sample based their 
flux at 8 and 24/im. 

When i nterpreted as an indication of evolutionary state the 77-2-25 index is at times know n to be 



misleading (JRobitaille et al. 



2006 



Crapsi et aU 2008; McCl ure et al. 



2010; 



Manoi et al. 



2011 



Edge 

on disks are known to have very red spectral indices (jPontoppidan et al.ll2005l ) and are frequently 
identified as Class I objects based on 712-25- The Spitzer IRS spectra reveal silicate emission and 
ice absorption features, which are not accounted for in an empirical 712-25 classification and allow 
us to gain a better understanding of the evolutionary state of these objects. Also important to 
note is K s band measurements from the extended 2MASS mission were taken in 2000 and the IRS 
spectra in 2007 and 2008. As these T Tauri stars are known to be highly variable this difference 
in epochs of observation can affect the 2 to 25/im slope and subsequent SED classification. The 
empirical SED classification from 2 to 25/j-m gives us one look into the evolutionary state of these 
objects. 
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4.2. Evolutionary State 

To better understand the evolutionary state of our sample, including those objects without 
known 2MASS K s band measurements, we examine the "ex tinction-free" spectral indices of the 
observed spectra for the entire sample (JMcClure et al.ll2010l ). These indices, 77,5-12 and ni2_20i 
characterize the SED slope of the spectra b efore extinction correction between wavelengths with 
roughly the same value of -r^ according to iMcClurd ( 20091 ). namely 5.3, 12.9 and 19.8/im. The 
flux at each of these wavelengths is calculated by averaging the flux over 3 channels, i.e. 0.3/um, 
centered at the anchor point. This index is truly extinction-free for objects with Ay>8 and varies 
slightly with wavelength for objects with Ay<8. Based on the 77,5-12 spectral index, we classify the 
evolutionary state of the sample (N = 79)according to the following criteria: 775-12 > —0.2 - the 
spectrum is dominated by envelope emission, —0.2 > 77.5-12 > —2.25 - dominated b y disk emission 
, and n5_i2 < —2.25 - dominated by photosphere emission iMcClure et al.1 ()2010l ). Using 715-12, 
23 objects are classified as envelope dominated, 55 as disk dominated, and 1 as a photosphere 
dominated. Large values of 77-12- 20, greater than 2.0, indicate typical values for TDs in Taurus and 
Ophiuchus (JMcClure et allboich . 



Figure plots these two extinction-free indices against each other, allowing us to gain another 
perspective of the evolutionary state of the objects in our sample. Three objects [#7, 35 and 137] 
fall within the area of this plotw here transitional disk s are found as they are characterized by a 



steeply rising SED beyond 12/mi (JMcClure et al.ll2010l ; iFurlan et al.ll201ll ). Upon examination of 



the SEDs of #7 and #35, the presence of deep silicate and CO2 ice absorption features in their 
spectrum an absence of detection in the three bands of 2MASS indicates these objects are deeply 
embedded in an envelope. Their decreasing SED from 5 to 12/xm and ste eply rising SED longer 
than ~10/iin suggests these deeply embedded protostars (JFurlan et al.ll2008l ) and are not TDs. The 
extinction-free indices are misleading here as the steep excess at longer wavelengths is a result of 
an envelope rat he r than, as would be the case for TDs, the excess originating at the outer disk wall 
(JEspaillatl 120091 ) . IFurlan et al.l (120 111 ) found embedded protostars seen edge-on are also found in 
this area of the plot. However, examination of the SED of #137 shows no evidence of a significant 
envelope affecting the spectrum, therefore the steep rise picked out in a high 77-12-20 is likely from 
emission originating in the outer layers of a heated outer disk wall. We will return to #137 in our 
discussion of radial structure in 



We examine the SEDs of our sample with respect to their 775-12 classification and identify 
objects which require more careful consideration and possible reclassification. Figure [7] displays the 
SED classes (based on 772-25) with respect to the extinction- free index 7x5-12- Of the 23 envelope 
dominated objects, 12 are Class I, 3 are FS and 4 are Class II. The other 4 lack data either 
at 25/i-m, as in the case of object #67, or in the K s band and therefore could not be classified. 
Ice and silicate absorption features, as a result of local extinction , are evidence for a s ignificant 



envelope affecting the SED (JBoogert et al.l 120081 : IFurlan et al. 



2008 



Zasowski et al.ll2009l ). We use 



the presence or absence of these features to get an idea of how strongly an envelope dominates the 
emission of these objects as well as whether their SED might truly be considered as dominated by 
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disk emission. The four Class II objects which are envelope dominated according to ri5_i2 [#68, 
69, 82 and 114] plus object #67 which has no ^2-25 but is envelope dominated based on 715-12, 
do not show ice or silicate absorption features in their observed or dereddened spectra; for these 
reasons we reclassify these objects as disk dominated though we do not include them in our disk 
analysis in § [5j If an envelope were present we would derive a large Ay for these objects in § [3j 
The relatively low extinction values calculated for #68, 69, 82 and 114 {Ay = 5.4, 3.7, 5.7, and 3.0 
respectively) signal an envelope does not strongly affect these systems. These objects are also close 
to the disk/envelope boundary in 715-12, especially #67, 68, 69, and 82 which have values between 
-0.20 and -0.09. This 'boundary' is tentative and an approximation based on the properties and 
locations of objects in Taurus and Ophiuchus in this plot. The uncertainties of #68, 69, and 82 span 
the disk/envelope boundary. Five Class I envelope dominated objects [#16, 19, 20, 113 and 115] 
display significant 10 /iin silicate emission features in their dereddened spectrum. Silicate emission 
features can be seen from the dusty disk in a Class II object, from the warm envelope in a face-on 
Class I object (jWatson et al.ll2004l ) and from both a disk and a warm residual envelope in a FS 
object. These objects appear more evolved than the canonical Class I stage and, though they were 
not classified observationally by n-2-25 as FS, these objects resemble the historical definition of a FS 
object as one which is in a transition phase between Clas s I and Class II, whose envelope is present 
yet t enuous and disk is beginning to be seen in the SED (jCasali fc Matthewslll992l : iGreene fc Lada 
19961 ). This explanation is consistent with the envelope dominated classification using n^-yi- 



Worth y of note in the enve lope dominated sample is object #122, or SSTc2d J032924. 1+311958, 
classified in iMerin et al.l (j2010l ) as a likely edge-on disk due to ice and slicate absorption features in 
the IRS spectrum. Our analysis classifies #122 as a Class I envelope. As stated in § 14.11 the 712-25 
classification of edge-on disks is known to be misleading. Ho wever the extinction free spectral index 
n 5-i2 classifies 6 of 7 edge-on disks in the Taurus sample of iFurlan et al.1 (|201ll ) as disks, with the 
one outlier being HARO 6-5B with an n^-yi = 0.26. Object #122 has an ri5_i2 of 0.69 as well 
as ice and slicate absorption features which suggest a significant envelope is still present in this 
system. Further, t he spectrum also resembles that of the famous edge-on disk, HH 30 in Taurus 
( Furlanet aDboilh . 



Of the 55 disks, 5 are Class I, 11 are FS and 36 are Class II according to 712-25- Three disk 
dominated objects [#3, 7, and 35] have no SED class as they were not detected in the K s band of 
2MASS ; as mentioned before, this is a strong sign these objects are deeply embedded protostars. 
The SEDs of #3, 7, and 35 display strong 10/jm silicate and 15/jm CO2 ice absorption features 
as well as a steeply rising SED beyond ~ 10/im. For these reasons we reclassify #3, 7, and 35 as 
envelope dominated. Three more objects classified as dominated by disk emission require closer 
examination [#13, 33, 36, and 63]. Class I based on 712-25, #63 has strong ice absorption features, 
a steep drop in flux from 8-9/im and a flat SED > 12/mi which is remin iscent of the most deeply 
embedded {Ay > 30) Class II objects in Ophiuchus (JMcClure et al.ll2010l ). The other objects, #13, 
33 and 36, are Class I and show strong silicate absorption and CO2 ice absorption features. These 
objects also seem to suffer from high extinction and demonstrate how the extinction-free indices 
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can give you a more accurate idea of evolutionary state, especially for objects with large Ay. Based 
on a closer examination of the n^-12 classification of objects in our sample with respect to their 
SEDs, our YSO sample in NGC 1333 consists of 21 objects dominated by envelope emission, 57 by 
disk emission, and 1 by photospheric emission; the SEDs of these objects are shown in Figure 
Figure [U and Figure [TUl respectively. 



According to lAspin et al.l (|1994l ). the photosphere dominated object #135 (ASR 54) has K = 
15.74 for which we would assign it to Class I. The SED of ASR 54 (Figure [TO]) shows 15/xm CO2 
ice absorption as well as ice absorption around 10/um, likely a result of f oregrou nd absorption and 
not the presence of an envelope. ASR 54 was noted bv lGutermuth et al.l (|2008bl ) as a candidate for 
a star which appears to be exciting PAH emission near the I RAS 4 filam e nt in t he southern NGC 



1333. F rom IRAC photometry and the K band detection by lAspin et al.l (11994! ) , iGutermuth et al. 
(|2008bl ) suggest this is a heavily extinguished (Ak ~ 10) late B type star with some excess at 



24/xm, impl ying it might be a transit ional or debris disk. The extinction-free spectral index ri5_i2 
agrees with IGutermuth et al.1 (j2008bl ) categorizing ASR 54 as a photosphere. There is certainly an 
excess at longer wavelengths, suggestive of a circumstellar accretion disk though we lack reliable 
extinction information which would allow us to determine whether this object is a transitional disk 
or debris disk. We also cannot conclusively say whether this is the object exciting the PAH emission 
seen in the southern end of NGC 1333. 



4.3. Class Fractions 



From our sample of full IRS spectra and 2MASS K s band photometry, the empirical 712-25 
Class I fraction is 23%±12% (16/70) and the fraction of sourc es dominated by envelope emission, 
as defined in the n 5 _i 2 analysis in §S2J is 27%±11% (21/79). IGutermuth et al.1 (ho08bh found 39 
of the 137 members of NGC 13 33 were protostars resulti ng in a protostar class fraction of 28%. 
Our findings are consistent with IGutermuth et al.l (|2008bl ) and reaffirms the protostar fraction is 
quite high in NGC 1333 which implies the presence of a substantial young population. 

In comparison to a cluster of similar age analyzed in the same way, such a s Ophiuchus (~0.8 
Myr, envelope fraction 9% (10/106), Class I fraction 25%: iMcClure et al.1 l201~oh . NGC 1333 has a 
significantly larger envelope fraction and a similar Class I fraction. Using Mathis's extinction c urve 
with Ry = 5, i.e. -^- = 3.06, the median Ay for Ophiuchus is 10.7 mag (IMcClure et al.ll2010l ) and 
the median Ay for NGC 1333 is 5.8 mag (see §[3j). While Ophiuchus is heavily extinguished, which 
skews age estimates from clas s fractions to the young er side, as extinction can make a disk have 
the SED class of an envelope (jPontoppidan et al.ll2005l ). NGC 1333 is not as affected by extinction 
and therefore its large envelope fraction demonstrates NGC 1333 is genuinely young. Given this, 
NGC 1333's larger envelo pe fraction would imp ly it is even younger than Ophiuchus whose age is 
estimated to be ~0.8 Mvr (|McClure et al.ll201( 



The fraction of our sample with disk dominated emission, as defined in the n§_i2 analysis in 
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§1121 is 72%±11% (57/79). Our sample of 70 objects with an SEP class yields a FS fr action of 
20%±12% (14/70) and a Class II fraction of 57%±12% (40/70). iGutermuth et al.1 (J2008bh identify 
98 Class II objects in their 137 YSO samp le, yielding a Class II fraction of 72%. Our results 
are consistent with IGutermuth et al.l (|2008bl )'s findings as their Class II sample is comparable to 
our disk dominated sample rather than our Class II sample as they do not classify objects as FS. 
Comparison with Ophiuchus, which has a larger fraction of objects with disk dominated emission 
(85%; 90/106), further demonstrates NGC 1333 is likely to be younger than Ophiuchus. 

Our selection for objects with excess in the 8/im IRAC and 24/mi MIPS bands does little to bias 
these class fractions as disk and envelope dominated objects typically have substan tial IR excesses 
at 8 and 24/im. In the study of NGC 1333 X-ra y sources by I Winston et al.l (|2010i ). they found 41 
Class III objects in their 95 object sample (43%). I Winston et al.l (120101 ) selected for likelihood of X- 
ray activity and therefore SED cl ass fractions would be bia sed toward later stages of star formation 
which have higher X-ray activity. IGutermuth et al.l (|2008bl ) calculated a disk fraction of NGC 1333 
using sources within a circle of radius of 5.5' or 0.4 pc centered at the median right ascension and 
declination of the YSOs they identified. Those authors calculate a disk fraction of 83%±11% as 
the fraction of the sample which have K s < 14 to all objects detected at K s after accounting for 
field star contamination. This method helps to mitigate biases against diskless members, which 
would make up the re mainder of the c luster specifically 17%±11%. In a study of low-mass stars 
and substellar objects, IWilking et al.l (J2004I ) estimate a disk fraction for the northern cluster of 
NGC 1333 to be 75%±20%, implying a diskless population of 25%±20%. They also use if-band 
excess for this calculation and correct for their incomplete census of members with disks. 



5. Disk Analysis 

In our analysis, we focus on the radial and vertical structure of the disk as well as the degree 
of processing of dust in the disks. We then statistically compare the structure of disks in NGC 
1333 to Taurus, Chamaeleon and Ophiuchus. The sample for our disk analysis consists of all disk 
dominated objects, according to ns_i2, for which we have properly characterized the extinction 
and which show no evidence of an envelope; this amounts to 44 of the 57 objects identified as disk 
dominated (8 FS and 36 Class II). 



5.1. Vertical Structure of Dust 



When a protostar forms it is surrounded by an infalling cloud of gas and dust which eventually 
lands as a flared accretion disk around the star. As the disk evolves, small dust grains coagulate, 
become larger, and eventually settle out of the upper layers of the disk down to the midplane. This 
decreases the amount of dust in the upper layers of the disk and reduces the heating and asso ciated 
flaring which causes the disk to become flatter (JD'Alessio et al.ll2006l ; iDullemond et al.ll2007l ). 
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The continuum emission from circumstellar accretion disks can be characterized through 
the use of spectral indices which sample silicate-feature-free emission, such as n.6-13 and ni3_3i 
(JD'Alessio et alJbood : iFurlan et allbood : IWatson et allbood : IFurlan et all 12009 ). The majority of 
the emission characterized by n6-i3 and rai3_3i is shown by iD'Alessio et al.1 ()2006l ) to originate 
in two distinct parts of the disk. Emission from the inner parts of a radially continuous disk is 
characterized by rt6_i3, which generally corresponds to 0.1-3 AU for the least settled disks and 
0.1-1 AU for the most settled disks. The middle to outer regions of radially continuous disks can 
be ch aracterized by n- 1 3-31 , around 3-25 AU in the least settled disks 1-10 AU for the most settled 
disks (JD'Alessio et alJbood ). 



We calculate n.6-13 an d ni3_3i for the our sample of disks in NGC 1333, as described above, 
using Equation [2] with fluxes at 6, 13, and 31^m derived by integrating the extinction corrected 
IRS spectra from 5.4-6.0/um, 12.8-14.0/iin, and 30.3-32.0/um respectively. The uncertainties in these 
spectral indices are propagated from the uncertainty in the original spectra. Some objects in 
the sample have significant uncertainty in their spectrum redward of 30/im resulting in a large 
uncertainty in ni3_3i. The increase in noise seen in some spectra at wavelengths longer than 30/um 
can be attributed to the data being taken late in the Spitzer mission (Campaigns 44, 48 and 49) 
after the detector bias was decreased in the LL array. The degree of settling in the disk can be 
inferred through comparison of these continuum indices with disk models which take into account 
settling. 

We compare our sample with predictions from models of irradiated accretion disks following 



the formalism of D'Alessio et al. 



f 1998. 119991 . 



2001 



2005 



20061 ). In the models used in this work, 



the surface density (and hence mass) of the disk is proportional to the mass accretion rate, M, over 
the Shakura-Sunyaev viscosity parameter, a. We chose to fix a at a typically adopted value of 0.01 
and vary the mass accretion rate. The reverse could have also been performed. Therefore, varying 
a is equivalent to varying the mass accretion rate. These models probe the effects of turbulence on 
the heating of the disk and thus the shape o f the SEP. Turbul e nt mixing has a non negligible effect 
over the vertical structure of the disk (i.e., IZsom et al.1 (1201 ll ). ICieslal (j2010l )). However, studying 
the role of the turbulence as a mixing mechanism and its effects on the emission of the disk is 
beyond the scope of this paper . Following the se models, a grid of disks around stars with M* = 0.2 
and 0.5 M Q was constructed (jEspaillati 120091 ) . The disks have a = 0.01, mass accretion rate M 
ranging from 10 -10 to 10~ 7 M Q yr _1 , and a dust depletion parameter e of 0.001, 0.01, 0.1 and 1 
(no settling). The continuum indices H6-13 and ni3_3i were computed for disk inclinations of 20°, 
40°, 60° and 80°. 

We compare these models with our observations in Fi gure [TT1 Also in Figure [TT] are the 
observed data f or Taurus (m i ddle panel; IFurlan et al.1 120061) and ou r disk sample in NGC 1333 
(bottom panel). IFurlan et al.1 ()2005l . 120061 ) and ID'Alessio et al.1 (120061 ) demonstrate as the settling 
parameter e decreases, n.6-13 and ni3_3i also decrease for the models described above, therefore 
lower values of n6-i3 and 7113-31 represent more settled disks overall. As discussed previously, 
it is important to remember these indices characterize the settling of dust in the upper layers of 
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different spatial regions of the disk. Overall the objects in NGC 1333 appear to be well explained 
by models with the e = 0.01 to 0.001. Dust depletion factors of 0.01 or 0.001 means 99 to 99.9% of 
the small grains, presumably from the ISM, have grown and sunk to the disk midplane, where they 
cannot be e asily seen. Th i s leads to an opt i cally ef fectively flatter disk, w ith a bluer spectrum, as 
discussed bv lFurlan et"afl (J2005I . bood . bood boilh : JD^lessio et~aD <|20od >. We find no statistically 
significant difference in the distribution of these continuum indices, as discussed in § 15.41 



5.2. Radial Structure 



5.2.1. Continuum indices and disk geometry 



Another continuum index, n 



, characterizes the excess emission from the inner most region of 



the disk (< 1 AU; lMcClure et al.ll2010i ). A disk which displays a large excess at longer wavelengths 
but has a very negative n<i-§ is similar to the SED of CoKu Tau/4 in Taurus, which has been 
modeled to have the geo metry of an inwardly truncated disk wit h a completely cleared out inner 
hole with radius ~10 AU ([D'Alessio et al.ll2005l ; ICalvet et al.ll2005l ). TDs, as defined in this analysis, 
have little to no excess from a disk at short wavelengths (< 8/nm) and show IR excess emission 
from an optically thick disk at longer wavelengths. PTDs, like UX Tau A, have an inner disk at 
<1 AU which would result in an short wavelength excess similar to that of a full disk though they 
posses one or more gaps and an outer optically thick disk. It is important to note that n<i-% does 
not necessarily characterize the amount of excess above the photosphere. 

The ni3_3i continuum index, as discuss ed previously, probes p roperties of dust in a radially 
continuous disk about 1-10 AU from the star ID'Alessio et al.1 (120061 ). For objects with interrupted 
radial structure, such as TDs, rai3_3i characterizes the excess originating from the outer disk wall, 
as the SED will show a steep rise at longer wavelengths which is essentially the superposition of 



a blackbody at the temperature of the d isk wall onto the spectrum of the star (JD'Alessio et al. 



2006 



Espaillat et al. 



2007 



Espaillatll2009l ). For TDs whose outer disk wall is heated due to stellar 



irradiation and accretion shocks on the stellar surf ace, we expect this warmer optically thic k dust 



to manifest itself in the spectrum as a high ni3_3i (ID'Alessio et al. 



2006 



Fur Ian et al. 



200J). 



I n this way, we can i dentify objects which have low ni— % and high ni3_3i as TDs (JMcClure et al. 



2010 



Manoj et al. 



2011 



PTDs in this analysis are disks with one or more gaps, resulting in at 
least an inner disk and an outer disk. Like TDs, PTDs with large gaps will have an increase in the 
continuum from 13-31/um from the wall of the outer optically thick disk. For this reason, a large 
^13-31 can identify PTDs with large (10s of AU) gaps, though PTDs with one or more smaller 



gaps would appear to have an ni3-3i similar to a rad ially continuous disk (JEspaillat et al. 



Espaillat 



2009 



Fur Ian et al. 



2009 



McClure et alJl20ld ). 



2007 



In Figure [TJl we plot ni3_3i vs. ni-% for the sample of disks as done in lMcClure et al.1 (J2010i ) 
to separate out objects whose radial structure is broken up by holes or gaps. As in Ophiuchus, 
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most of the disks in NGC 1333 fall within -2 < n2-6 < -0.7 and -1.33 < ^13-31 < 0.3. In Ophiuchus 
and Taurus, PTDs and TDs are found in the upper octile in m 3-31 with TDs in the lower octile 
of ri2-6 an d PTDs between the lower octile and the median (JMcClure et al.ll2010l ). Applying this 
classification to NGC 1333, #137, the outlier in the TD range of ni2-20> a l so stands out here in 
^13-31 an d has a very steep slope from 2-6/jm placing it in the TD range as defined by Ophiuchus 
and Taurus of the ni3_3i vs. ri2-§ plot. Six objects [#55, 75, 101, 110, 131, and 134] fall within or 
very close to the supposed PTD range of ri2_6, as they are outliers in ni3_3i- All of these objects 
appear to be quite faint from 1-35/xm with v¥ u ~ 10" 11 to 10~ 12 ergs s _1 cm -2 . Two objects, #55 
and 75, do not appear to have a rising SED at longer wavelengths but instead have a ni3_3i outlier 
as a result of noise in the spectrum > 30/um. The three objects in the PTD range with the largest 
values of ni3_3i, #101, 110, and 131, appear to have SEDs similar to TDs, as they show steep 
slopes from 2-6^m and flux excess from 13-31/xm. However, #131 has an excess <8^m similar to 
that of a disk and is likely not a TD. These objects do have a steep slope from 2-6/xm, with U2-q 
~ -1.7. A small amount of dust in a central clearing would cause the slope to be slightly flatter, 
steeper than a full disk but flatter than a photosphere. For these reasons we consider #101 and 
110 as possible TDs. 

The obje ct with the lowest value of Wjs-aj in the PTD range, #134, does not look like any 
known PTDs (|Espaillatll2009l ; lEspaillat et al.ll2010l ) but it shows a short wavelength excess similar 
to a disk, a significant dip in flux between 10 and 20/im and starts to rise out to about 30/um. 
We do not classify #134 as a TD or PTD though note the strong dip in the middle of the IRS 
spectrum possibly indicates interrupted radial structure and makes this object a good candidate 
for further study. Iwinston et all (J20ld ) identified #136 as a TD based on its IRAC and MIPS 
photometry. Though it is not an outlier in ni3_3i, #136 has the steepest 712-6 (-2.6) which supports 
its classification as a TD. We will discuss #136 further in § 15.2.21 As mentioned previously n2_6 
is a slope and does not necessarily characterize the amount excess above the photosphere. For this 
reason it may not be able to identify some types of PTDs and TDs. 

There are two blue outliers [#53 and 123] in ni3_3i with an rii3_3i below -4/3, the limit 
for an infinite geometrically thin, opti cally thick disk. Low values of ni3_3i have been found in 
models of outwardly truncated disks (JMcClure et al.l 120081 ) and for this reason we suggest #53 
and #123 have outwardly truncated disks. Important to note is objects may also have an ni3_3i 
below -4/3 if they are no longer optically thick, a sign of a more evolved state. However, we do 
not believe this is the case for these objects as they are not identified as Class III from § 14.11 or 
photosphere dominated from § 14.21 The presence of a companion i s known to truncate a disk (e.g. 
Artvmowicz &, Lubowlll994l ) as in the case of SR20 in Ophiuchus (JMcClure et al.ll2008l ). The lack 
of established multiplicity information for the sample does not allow us to address definitively the 
multiplicity for these truncated disks in NGC 1333. 
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5.2.2. Silicate emission and disk geometry 

The broad silicate emission features in the spectrum of disks at 10 and 20/im arise from 
optically thin dus t at the disk surface, and can be used to probe the dust pro perties in the upper 



layers of the disk (JMalfait et al. 



1998 



Natta et al. 



2000 



D'Alessio et al.ll2006l ). For objects whose 



structure deviates from being radially continuous, such as TDs, the emission seen in these features 
comes from optically thin dust within the gap or hole, such as in GM Aur, or from optically thi n 
dust in the surface layers of the outer disk wall, such as in CoKu Tau/4 (JD'Alessio et al.l 120051 ) . 
Therefore measuring the strength of the 10 and 20^m silicate emission can also be used to probe 
the radial structure of the disk. 

The equivalent width of the silicate emission features in the IRS spectra is a measure of the 
optically thin emission per unit area of optically thick emission from dust in the disk. The equivalent 
width is characterized by the sum of the excess flux from the silicate feature at each wavelength, 
F\, with respect to the flux of the underlying continuum, F\^ cont . Equivalent width is given by 



EW X 



A 2 



■r\ r \,cont 
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X,cont 



(3) 



In our analysis we adopt Ai = 8^m and A2 = 13/xm and Ai = 16/im and A2 = 28^m in order to 
measure the 10 and 20//m silicate emission features respectively. We adopt positive values of EWa 
for emission features. 

While the equivalent width expresses the strength of the silicate feature in units of its under- 
lying continuum, the integrated continuum subtracted flux of the silicate feature allows us a probe 
how much optically thin dust is actually producing this emission. For this reason, we also calculate 
the flux of the 10 and 20//m silicate features using 



F> 



A 2 



Ai 



(F\ — F\^ cont )d\ 



(4) 



with the sam e limit s as above. The EW and integrated flux a r e calculated i n the same way as 
Watson et all I^POsI ): iFurlan et al.l |200g| ): JMcClure et al.l |201(il ): iManoi et"aD l|201ll ). 



We determine the continuum for the calculation of both properties by fitting the spectrum with 
a 3rd to 5th order polynomial with the fit anchored at 5 . 61-7.9 4//m, 13. 02-13. 50/um, 14. 32-14. 80/iin, 
30. 16-32. 19/im and 35. 07-35. 92^m which iWatson et al.l (J2.009) showed to be free of distinct silicate 
emission features. For objects where a single polynomial did not produce an adequate fit, we used 
two polynomials to fit the continuum from 5-14//m and 14-36/im. For a few cases a smoothly 
varying fit which did not exceed the extinction corrected flux of the spectrum could not be found. 
For these objects the equivalent width and integrated flux from the portion of the spectrum for 
which a satisfactory fit could not found are left out of our analysis. We calculate EW10, EW20J E10 
and F20 as well as their uncertainties for 44 of the 57 disk dominated objects with no evidence of an 
envelope as discussed in §[57TJ these values are listed in Table 01 The main source of uncertainty in 
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these calculations results from the continuum fit, the uncertainties noted in Table [Hare calculated 
assuming an uncertainty in the continuum fit of 10%, though in some cases the uncertainty is likely 
to be larger. A continuum fit successfully extracting the 33/im crystalline forsterite feature could 
not be found for all but a few of the spectra as a result of significant uncertainty at wavelengths > 
30/im as discussed in § 15.11 therefore we do not discuss this feature in our analysis. 



The EW c an be used to iden tify PTDs which have gaps which result in a deficit of optically 
thick emission (JFurlan et al.ll2009l ). Equivalent width measures optically thin emission with respect 
to optically thick emission, therefore for the same amount of optically thin emission, a lack of 
optically thick emission would manifest itself in a large EW. If the gaps have some amount of 
optically thin dust this will further enhance their ability to be identified by h aving a large EW. I n 
Figure [13] we plot ni3_3i vs. EWio for our sample and the Taurus sample from lFurlan et al.l 1)20061 ) . 
along with a polygon indicating the boundaries in these indices which correspond to boundaries 
found by the models of radially continuous, irradiated disks described in S I5.11 These boundaries are 
a function of the parameters varied in the models (e.g., M, M and disk inclination). We also note 
that the EWio depends on the settling parameter^ which we also varied in the models. Objects 
with EWio and ni3_3i which fall outside of this polygon have an interrupted radial structure, 
meaning the disk has one or more gaps or is inwardly or outwardly truncated. This is reflected the 
positions of the well modeled TDs and PTDs in Taurus such as CoKu Tau/4, GM Aur and LkCa 
15. Figure [HI demonstrates the areas outside the polygon of radially continuous disk models which 
generally correspond to dis ks with altered radial structure in Taurus, Chamaeleon I an d Ophiuchus 
which have been modeled (JEspaillat et al.l I2OO7I ; HVlcClure et all hood ; kirn et al]l200a ). 



First, we examine the disks which are extreme outliers in EWio through the criterion that 
their EWio is in the top 12.5%. Using this restrictive criterion we identify #50, 52, 57, 73 and 
116 as PTDs. As discussed previously, objects which are likely to be TDs possess a large ni3_3i, 
greater than that of the radially continuous disk models. While #110 falls within the PTD range 
in Figure [12] it is certainly not an outlier in EWio, though it is an outlier in ni3_3i reinforcing its 
previous classification as a TD. Curiously, #110 has a position in Figure [T3~l similar to UX Tau A, 
a PTD which has bee n modeled to have no optically thin dust in the gap resulting in a small EWio 
(JEspaillat et al.ll2007l ). We cannot confirm nor exclude this possibility without detailed modeling. 
All disks which were in the PTD range of Figure [T2l except for #110, fall within the polygon of 
the models in Figure [T3l The other TD in the PTD range from Figure [T2l #101, falls just within 
the polygon of radially continuous, irradiated disk models in Figure [13] though its uncertainty in 
n i3-3i crosses the top boundary. Of the other disks in the PTD range in Figure [T2l #75 and 134 
have EWio greater than the median, though are not in the upper 12.5%, therefore are not large 
enough to be considered significant outliers. The classification of #137 as a TD from Figure [6] and 
Figure [12] is bolstered here as its position in Figure [13] places it outside the polygon of radially 
continuous disk models and in the TD range. Our analysis agrees with that of lMerin et al.l (]2010l ) 
which classifies #137 (known therein as source 5 or SSTc2d J032929. 3+311835) as a 'cold disk.' 



The TD suggested by IWinston et al.l (]20ipJ) which did not fall in the TD range of Figure [T2] 



20 



is also not unique here in it falls within the polygon of models. Though this object has a 'normal' 
ni3_3i, examination of it SED shows it is a TD. The SED of #136 is photospheric < 9/jm and the 
continuum slightly rises toward the long wavelength side of the spectrum. We interpret this as a 
sign of a large central hole size, relatively cool outer disk wall, and a settled low mass outer disk. 
Both the large inner hole and the low-mass settled outer disk prevent #136 from being identified 
as a TD based on ni3_3i- 

The case of #136 suggests ni3_3i can pick out some TDs though others can have an ni3_3i 
similar to that of a radially continuous disk. The spectral indices ny^^i and n<i-§ are only good at 
identifying TDs with a certain range of hole sizes. If the hole is too small (<1 AU) or too big (>50 
AU) the indices will not identify them as TDs. It is also possible ^13-31 can improperly identify 
radially continuous disks as TDs therefore individual SED inspection is key for proper classification. 
Further discussion of the classification TDs and PTDs in NGC 1333 is found in § [6j 



5.2.3. Analysis and comparison of silicate emission features 

With respect to this classification we further analyze the radial structure of disks in NGC 1333 
using Fio, EWio, and EW20- Figure [TBI shows Fi (top) and EW20 (bottom) plotted against EWio- 
As expected, Fio and EWio are we U correlated, with Pearson's linear correlation coefficient r = 
0.65, and a probability p vn ,nd = 0.0024% that this correlation could have been generated from a 



random distribution (see IWatson et al.l (|2009l )). A small value of p ran d, namely < 1%, indicates 
the quantities are significantly correlated. The Fio vs - EWio plot shows most of the PTDs have 
large values of Fio as well as an enhanced EWio- A handful of disks have similar F 10 values as the 



PTD s but have much s maller EWiq. Obje cts with enhanced EWio i n Ophiuchus (IMcClure et al 



2010 ) and Chamaeleon I (|Manoj et al.ll201ll ) show similar values of Fio to radially continuous disks, 
which suggests these objects have a deficit in optically thick emission as a result of a gap in their 
disk. Interestingly, the disk with the highest Fio also has the highest EWio i n NGC 1333, which 
is not seen in Ophiuchus and Chamaeleon I. The TDs in NGC 1333 have low Fio, also the case for 
TDs in Ophiuchus and Chamaeleon I. The flux at 10/im comes from the inner regions of the disk, 
therefore a significantly cleared inner disk, as is the case for TDs, would result in a lack of emission 
at 10/iin. 

The bottom panel of Figure HU compares EWio and EW20 and shows the two trace each other 
well, with r = 0.32 and p r and = 6.1%. All of the PTDs, which have enhanced EWio, fall above the 
median for EW20- However most of the PTDs have a lower EW20 than EWio which is also seen 
in Chamaeleon I. A large gap would affect the underlying optically thick emission contributing to 
EW20 as well as EWio- This trend in lower EW20 with EWio outliers possibly suggests most of 
these EWio outliers have small gaps in the inner disk which only affect the EWio- As discussed in 
§ I5.2.H the small nature of these gaps would cause them to be more difficult to identify as PTDs 
with the n2_6 vs. ni3_3i analysis (Figure fl"2l) . As was the case in Ophiuchus and Chamaeleon I, 
the trends of the continuum indices, fluxes and flux ratios with EW20 are weaker versions of the 
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trends with EWio- Most of the scatter in EW20 is due to the uncertainty in the continuum fit for 
the 20/mi feature. 



5.3. Dust Processing 

One probe of the degree of crystallization and grain growth in the inner 1-2 AU of the disk is the 
flux density ratio Fi 1.3 /Fog (IPrzygodda et all 20031; Ivan Boekel et al.112003 . 12003 : iKessler-Silacci et al 
20061 : iHonda et al.ll2006l ; iBouwman et al.ll2008l ; IOlofsson et al.ll2009l ). F11.3/F9.8 quantifies the shape 
of the 10/im silicate feature and as such tracks the degree of processing of the dust which creates the 
silicate emission features. The higher the ratio, the more the dust has been processed (more grain 
growth and/or crystallization), with a 'pristine' profile defined as one similar to the optically-thin 
interstellar-grain profile with F11.3/F9.8 = 0.34. We calculate Fn.3 and Fg.g according to Equation 
4 with Ai = 10.9 and 9.4/iin and A2 = 11.7 and 10.2/mi respectively. F11.3/F9.8 for the sample is 
listed in Tabled! An F11.3/F9.8 representative of the processing of dust in the disk of object #110 
was impossible to obtain due to contamination of Fg.g by residual molecular hydrogen emission not 
removed by sky subtraction. 

In Figure[16]we plot F11.3/F9.8 against EWio- The objects with the highest values of F11.3/F9.8 
correspond to the lowest EW strength, which is characteristic of radially continuous disks. This 
inverse relationship has Pearson r = -0.34 and p ran d = 2.4%. We also notice the trend that as 
EWm decreases the frequency of disks with high Fn.3 /F9.8 increases. This trend is seen in Ta urus 



( Furlan et al.l 12000 ) . Qp hiuchus (McClure et al.1 201C) and Chamaeleon I ( Manoi et al. 



well as in other regions (Ivan B oekel et al.l |2005|; IKessler-Silacci et al.l 120061 : 



2011 



Bouwman et al. 



), 



as 



2008; 



Apai et alJl200a ; IPascucci et al.ll2009i ). Also, the PTDs have relatively low F11.3/F9.8, with most 



of them falling below the median and close to the pristine line. This would imply in general the 
dust seen in the PTDs has experienced less grain growth and/or crystallization than the radially 
continuous disks. Therefore the optically thin dust at 1-2 AU in PTDs is more pristine than 
radially continuous disks. The trend in low F n.3 /Fg.g (more pristine shape) in outli ers in EWm is 



also seen in model s of PTDs (lEspaillatl 12009) and observati on of PTDs in T aurus (jSargent et al. 
20091 ). Ophiuchus (JMcClure et alJboioh and Chamaeleon I (JKim et alihood ). Modeling suggests 



the emission responsible for the shape and strength of the 10/im silicate feature in these di sks comes 
from the optically thin dust within the gaps or in the surface layers of the outer disk wall (jEspaillat 
2009). As also seen in Chamaeleon I, the disks with ni3_3i below -4/3 in NGC 1333 show a slight 
preference toward higher values of F11.3/F9.8 which implies these truncated or optically thin disks 
have a higher degree of dust processing than TDs or PTDs. 
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5.4. Statistical comparison to Taurus, Ophiuchus and Chamaeleon I 

We use a one dimensional two-sided Kolmogorov-Smirnov (K-S) test to determine if the distri- 
bution of n2-6, rc-6-1 3; and rei3_3i in N GC 1333 is draw n from a significan tly different population 
from that o f Taurus (IFurlan et al.ll2006l ). Chamaeleon I ( Manoj et al.ll201ll ). or the L1688 cloud in 



Ophiuchus (IMcClure et al.ll201Cll ). To do this we took our sample of 36 Class II disk dominated 
objects and did not include objects which were identified previously as TDs, truncated disks or 
PTDs, leaving us with 27 objects representative of the 'normal' radially continuous disk population 
in NGC 1333. We culled the data from Taurus, Chamaeleon I and L1688 in the same way. The 
distribution of these indices is show in Figure H7] and p- values inset from a ID two-sided K-S test 
comparing NGC 1333 to each region. 

We find no statistically significant difference between NGC 1333 and Chamaeleon I, Taurus, 
or L1688 in the distribution of ng_i3 and 7113—31. These statistical tests imply there is a similar 
dispersion in these continuum indices in all four regions. The distribution implies a heavy concen- 
tration in the e = 0.001-0.01 range and a tail that reaches 0.1, but very few objects that demand 
a well mixed disk. We do, also, find a difference between NGC 1333 and Taurus and L1688 in 
the distribution of 712-6 (p-value = 0.07 and 0.01, respectively). This would imply there are sig- 
nificant differences in the slope of the excess emission in the innermost regions of disks in these 
regions. NGC 1333 and Chamaeleon I have distributions of U2-6 skewed toward bluer (smaller) 
values. In Chamae leon I, this is attributed to a spectral type difference between these two regions 
( Manoi et al.ll201ll ). The same argument can be made for NGC 1333 which has many more objects 
in the M3-M8 than K5-M2 spectral type range and a median spectral type of M3. Taurus and 
Ophiuchus both have a median spectral type of M0. The median U2~q value for M3-M8 objects in 
NGC 1333 is -1.64 which when compared to the median ni-% for K5-M2 objects of -1.50 shows the 
prevalence of M3-M8 stars in or sample skews the ri2-6 distribution of NGC 1333 to the smaller 
side. 

Trends with model parameters within the polygon from Figure [T3l are shown with arrows in 
Figure H (JFurlan et all bood ; b'Alessio et al.1 bood : lEspaillatl bood ). To determine if there are 
difference in the scatter of points within this polygon a 2-dimensional two-sided K-S test was 
performed for only the disks in Taurus, Chamaeleon I, Ophiuchus, and NGC 1333 which fall within 
this polygon and have a spectral type later than M0. We find no statistically significant difference 
in the 2D distribution of points within this polygon between NGC 1333 and Taurus, Chamaeleon 
I, or L1688, with p-values of 0.36, 0.20, and 0.21 respectively. 



5.5. Median Spectra 

We construct a representative median for our sample in NGC 1333 from the 35 objects classified 
as both Class II and disk dominated observed in program 40525. The objects in our median span 
a spectral type range from A3 to M8, with 7 between K4 and M2 and 20 between M3 and M8 
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and a median spectral type of M4. We also include 7 objects w ith unknown s pectral type. The 
extinction measured in Ay (derived from the Ry = 5 version of iMathid (|1990l ) extinction curve) 
ranges from to 22 .3, with a me d ian of 3.1 and includes seven objects with Ay> 10. We follow the 
same procedure as iFurlan et al.l (|2009i ) to construct our median. We normalized the dereddened 
IRS spectrum of all 35 Class II disks through multiplication by a scalar factor equal to the median 
dereddened H-band flux of the 35 objects (0.039 Jy) divided by the i^-band flux of each object, 
computed the median of the sample at each wavelength and converted to v^ v . This method 
effectively normalizes the excess emission for differences in stellar luminosity for stars with spectral 
type K-M (JD'Alessio et al,lll999i ). which comprise the majority of our median sample (27/28 objects 
with known spectral type). 

In Figure [18] we compare our media n to the medians of K5 to M2 stars in Taurus, LI 688, 
and Chamaeleon I previously presented in lFurlan et al.l (J2009J) each scaled to the median .ff-band 
flux of Taurus. The shape of the NGC 1333 median is quite similar to that of the other medians, 
across the IRS range. After scaling to the median H-band flux of Taurus, the NGC 1333 median is 
still fainter than the others. This can be attributed to the spectral type differences between these 
samples. The Taurus, L1688, and Chamaeleon I medians are for objects with K5 to M2 spectral 
types while our sample has six objects in the K5 to M2 range and 20 have spectral types later 
than M2. Late type stars heat their disks less as they are less luminous and therefore irradiate 
their disks less as well as have lower mass accretion rate resulting in weaker accretion heating. The 
cooler nature of late type stars manifests itself as fainter IR excess emission. 

We compare the shape of the 10/iin silicate features of the medians of Taurus, L168 8, Chamaeleon 



I, and NGC 1333 in Figure QjJl We do this in the ex act same way as de scribed in IFurlan et al. 



(J20091 ) . The only difference between our Figure [19] and IFurlan et al.l (|2009l ) 's Figure 7 is the addi- 
tion of the 10/mi silicate feature of NGC 1333 to the plot. Also plotted in Figure [T9l is a 'pristine', 
ISM-like profile which is represented by the average of t he continuum-subtr acted and normalized 
10/iin silicate feature of LkCa 15 and GM Aur in Taurus (jSargent et al.ll2009i ). This pristine profile 
is scaled to the mean 9.8/^m flux of the four regions. All four regions show evidence for grain growth 
and/or crystallization as they show a silicate profile wider than the 'pristine' profile for ISM-like 
small amorphous silicate grains (sub-/im). NGC 1333 shows a long wavelength wing similar to that 
of the other regions which suggests similar grain sizes in the optically thin dust. 



6. Discussion of PTDs and TDs 



We have examined two ways of identifying TDs and PTDs from properties of their spectra in 
§ 15.2.11 and § 15.2.21 The first analysis, that of the ni3_3i vs. U2~q diagram (Figure [T2|) . identified 
outliers in rai3_3i, all but one of which fell within the PTD range as defined by the positions of 
TDs and PTDs in Taurus and Ophiuchus. Secondly, we analyzed the equivalent width of the lO/xm 
silicate emission feature and identified strong outliers. An enhanced EWio is suggestive of a deficit 
in optically thick emission relative to optically thin emission, as would be the case if there are gaps 
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in the disk. We also used the limits of model indices in the ni3_3i vs. EWio diagram (Figure [T 
to examine the PTDs and TDs from the first analysis. All of the outliers from the first analysis 
except #110 and #137 fell within the boundaries of model disks, though th e uncertainty in m.s 3 1 
for #101 does cross out of the range of model disks. The TD suggested bv lWinston et al.l (|2010l ). 
#136, did not stand out in any of the previous analyses, except it had the most negative n 2 -e 
which is closest to that of a photosphere. 

We now examine the TDs and PTDs in a third way, with respect to a representative median 
of the population of disks in NGC 1333 from § 15.51 In this analysis, we identify PTDs as having 
near-IR excess at or below the representative median of NGC 1333 scaled to the object's H-hand 
flux, strong silicate emission with flux exceeding that of the scaled median, and flux comparable 
to the scaled median longwards of 24^m. While object #73 exceeds the flux of the median at long 
wavelengths, it s spectrum resembles known PTDs such as ROX 44 (JEspaillatll2009l ). as well as the 
"kink" disks of lSicilia-Aguilar et al.l (|201ll ). for these reasons we classify it as a PTD. We identify 
TDs as having a flux deficit with respect to the scaled median in the near-IR and flux exceeding the 
median at wavelengths greater than 24/im. It is important to note that the dependence of these 
classifications on emission around 24^m restricts the range of hole s i zes p resent in transitional 
objects identified in this manner. The TD identified bv lWinston et al.l (|2010i ). #136, notably does 
not conform to this definition as it falls below the scaled median for the entire IRS range, however 
it is clear that this spectrum is representative of a disk transitioning from a the spectrum of a full 
optically thick disk to that of an optically thin debris disk, and therefore we classify it as a TD. 
Object #116 cannot be classified with respect to a scaled representative median of the disks in 
NGC 1333 as its flux is much greater than the scaled median. 

Figure [20] shows the dereddened SED of the 9 PTDs and TDs along with the representative 
median of NGC 1333 scaled to each object's -ff-band flux. Table [5] shows what classification each 
analysis suggested and our final classification of these disks. We favored analysis of the SED with 
respect to the median as it has been reliably used to identify TDs and PTDs in many previous 



studies (ID'Alessio et al 



Kim et al 



2009; 



2005 



Merin et al. 



Brown et al.1 120071 : ICalvet et al.1 1200.4 lEspaillat et al.1 120071 . 12003; 
20101 ) . Also, as stated § 15.2.11 and § 15.2.21 analyses using ?i2_6 and 



ra i3-3i are only sensitive to a certain range of hole or gap sizes and will not effectively pick out small 
holes or gaps. In most cases, classification based on the EWio vs. ni.3-31 agre e d wit h the SED 
analysis. In all cases the ni3_3i vs. ?i2_6 analysis, as defined by iMcClure et alJ (J2010l ) disagreed 
with the classification from the other two methods. 

At least 9 of the 44 disks analyzed in our sample of NGC 1333, 20%, show evidence for gaps or 
holes in their disks in various analyses. In our sample we identify 4 as TDs, which translates to 11% 
of disks analyzed in our sample NGC 1333 (4/44). We do not have sufficient multiplicity information 
for our sample to address the contributions of companions on our transitional disk fraction. NGC 
1333's TD fraction is larger than the transitional disk fractions in t he 1-2 Myr Taurus region (3.5%), 
~2 Myr Chamaeleon I region (5.8%) and L1688 ~0.8 Myr (3.2%) (JRirlan et al.ll2009l ). Though the 
sample of disks we analyzed is relatively small (N = 44) , observations indicating gaps in disks are 
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suggestive of giant p l anet formation reg ardless of sample size (JFurlan et al.ll2009l ; iKim et al J 12003 : 



McClure et al. 



2010; iManoi et al.ll201ll ). At the young age of <1 My r, disks in NGC 1333 are 



forming planets and dissipating their disks. Furthermore, lAspinl (|2003l ) cite age estimates for 4 of 
the PTDs identified in our analysis, #50, 52, 57, and 116, with ages of 1 Myr, 0.5 Myr, 3 Myr and 
1 Myr, respectively. It is important to remember that age estimates of young stars using isochrones 
are known to be uncertain as discussed in § [TJ However, the dispersion in these ages might suggest 
the PTD phase can occur more than once in the life of a disk. 

The finding of PTDs at very young ages supports theories for the rapid formation of giant 
planets such as through core accretion with osci llatory migration (IRice &; Armitagd 120 03) or grav- 



itational collapse of gas and dust in the disk (JGoldreich fc Ward 



1973 



Boss 



200l| ). The quick 



timescale for planet formation demonstrate d in NGC 1333 is not currentl y accounted for in some 
models of planet formation such as thos e in Alexand er fe; Armitagd (120091 ) or in the older core ac- 



cretion models of giant planet formation (jPollack et al. 



w- 



1996). However, our results do not discount 



theorie s of "born again" disks which suggest g iant planets could form more than once in the life of 



a disk (jSargent et al. 



2006 



Watson et al.l 120091 ) . 



7. Summary and Conclusions 

We presented the IRS spectra of 79 YSOs in NGC 1333. We supplemented the IRS spectra 
with published photometry and constructed SEDs. We used the SEDs, before extinction correc- 
tion, to classify our sample into SED classes as well as evolutionary states, using "extinction-free" 
indices. Based on this classification, our sample in NGC 1333 consists of 21 objects with envelope 
dominated emission, 57 objects with disk dominated emission and 1 object with photosphere dom- 
inated emission. Further, we analyzed various spectral properties of 44 objects which show disk 
dominated emission and had no evidence of an envelope in their observed SEDs. We compared the 
shape of the continuum emission as well as silicate emission features to models of radially continu- 
ous irradiated accretion disks in an effort to understand the vertical and radial structure of these 
disks in NGC 1333. We used two different methods in order to identify disks with altered radial 
structure. With these analyses, as well as comparison with a representative median spectrum, we 
classified TD and PTDs in our sample. We also compared the properties of disks in NGC 1333 
to disks in the Taurus, Ophiuchus, and Chamaeleon I star forming regions in order to investigate 
potential differences with age and population. Our main results are summarized below: 



We find no difference between the distribution of dust settling continuum indices, ri6-i3 
and 7ii3_3i in NGC 1333, Taurus, L1688 in Ophiuchus, and Chamaeleon I. However, through 
comparison of class fractions, we find NGC 1333 is likely to be even younger than the youngest 
of those three regions, Ophiuchus. This similarity between star forming regions of significantly 
different ages implies advanced dust settling (e ~ 0.001) can happen very early in the life of 
a disk. Alternatively, few disks are well mixed even at the young age of < 1 Myr. 
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• Of the disks analyzed in NGC 1333, 20% show evidence of radial gaps and clearings and 
11% are identified as transitional disks. These fractions are significantly larger than those of 
Taurus, Ophichus and Chamaeleon I. We find these observations of NGC 1333 to be evidence 
of giant planet formation resulting in disk dissipation at a very young age, < 1 Myr. 

• Most PTDs, which are large outliers in EWio, were not large outliers in EW20- This suggests 
the gaps in these disks which produce an enhanced EWio are n °t large enough to affect EW20- 
The smaller gaps in these PTDs likely results in an ni3_3i similar to radially continuous disks 
and makes them difficult to identify in analyses which use a high ni3_3i as a criterion. 

• Multiple analyses are needed to identify PTDs from spectra without detailed modeling. The 
ri2-6 vs. ni3_3i was not effective in selecting PTDs in NGC 1333 and in some cases con- 
tradicted the classification via EWio y s- ^13-31 and comparison to a scaled representative 
median. 

• Comparison of a representative median of disks in NGC 1333 to Taurus, Chamaeleon and 
L1688 demonstrates the median spectrum of NGC 1333 had a similar shape to all three 
regions, though after scaling was still fainter. The prevalence of late type stars in our NGC 
1333 sample can account for this flux discrepancy. 

• PTDs in our sample NGC 1333 have more pristine optically thin dust in the inner parts of 
their disks based on analysis of the shape of their 10/um emission using F11.3/F9.8. The dust 
in these disks with radial gaps seems to have experienced less grain growth and crystallization 
than radially continuous disks. 



This work is based on observations made with the Spitzer Space Telescope, which is operated 
by the Jet Propulsion Laboratory (JPL), California Institute of Technology (Caltech), under the 
National Aeronautics and Space Administration (NASA) contract 1407. In this work we make use 
of data products from the Two Micron All Sky Survey, which is a joint project of the University of 
Massachusetts and the Infrared Processing and Analysis Center/California Institute of Technology, 
funded by NASA and the National Science Foundation. This research also utilizes the Vizier and 
SIMBAD databases, operated at CDS, Strasbourg, France, NASA's Astrophysics Data System 
Abstract Service. Support for this work was provided by NASA through contract number 1257184 
issued by JPL/Caltech, JPL contract 960803 to Cornell University, and Cornell subcontracts 31419- 
5714 to the University of Rochester. 
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Fig. 1. — Spitzer IRAC, MIPS and IRS data on emission from the southern outflow lobe of IRAS 
4B. The S(7) to S(l) H2 lines are identified. The LL spectrum contains many attenuated blends of 
H2O emission lines from IRAS 4B, which li es 5" north of t his po sition. The IRS spectrum of IRAS 
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analysis. 
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Fig. 2. — SED of a background galaxy redshift of 0.29. The dotted lines indicate the 6.2, 7.7, 8.6 
and 11.2/xm PAH features at their redshifted wavelengths. 
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Fig. 8. — Continued... 
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Fig. 9. — SEDs of objects dominated by disk emission, as defined in the n.5-12 analysis in § 14.2 
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Fig. 10. — Observed SED of photospheric object aka ASR 54. Symbols and lines are the same as 
in Figure 
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Fig. 11. — rt6-i3 vs. ni3_3i indices for radially continuous model accretion disks ("see § 15.11 for 
details), Taurus (indices from P. Manoj, personal communication using spectra from iFurlan et al.l 
(2Q06| )) and NGC 1333. The triangle indicates where n@_i3 and ni3_3i are that of a infinite 
geometrically flat optically thick disk. The value for a photosphere for «6-i3 and ni3_3i is -3 and 
plotted with a star. Similar to Figure 22 from iMcClure et al.l (|2010l ). Typical uncertainties for 
Taurus and NGC 1333 are shown below the region's name. 
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Fig. 13. — EWio vs - n i3-3i for our sam ple of disks in Tau rus (top, indices from P. Manoj, 
personal communication using spectra from iFurlan et al.1 ( 20061 )) and NGC 1333 (bottom). TDs 
are plotted here with squares, PTDs with stars and truncated disks with diamonds. The dotted 
line indicates the median value and the dashed line indicates the upper octile. The polygon denotes 
the boundaries of these indices in models of radially continuous irradiated accretion disks (see § 15.11 
for details). 
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thick disk (-4/3, indicated with a dashed line). This trend demonstrates areas of the polygon which 
correspond to disks which are well mixed or settled. The light gray arrow shows how varying the 
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Fig. 15. — Distribution of Fio and EW20 y s- EW10 for objects in NGC 1333, symbols are the same 
as Figure [13j Typical uncertainties are shown in the upper right hand corner of each plot. The 
dotted lines indicate the median value. The dashed line in the bottom plot indicates where EW10 
and EW20 are equal. 
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Fig. 16. — Distribution of the flux density ratio F11.3/F9.8 vs. EW10 for objects in NGC 1333, 
symbols are the same as Figure [T3j The dotted line indicates the median value and the dashed line 
indicates 0.34, the value of F11.3/F9.8 of pristine dust. 



54 



30 
25 
20 

15 

10 

5 

3§ 

30 

w 

O 25 

" 20 

O 15 

O 10 

z 

5 
38 
25 
20 
15 
10 

5 



'2-6 



£ 




Taurus 
NGC 1333 



.07 :: 



"I I 1 ' 



'2-6 



Chamaeleon 
NGC 1333 




.92 :: 



'I 1 I" 



L1688 

" l 2-6 NGC 1333 



3 



i£A~ 



P = .01 



'6-13 



:: n 



6-13 



P = .96 



XI 



r. 



11 n 



6-13 



P = .98 






c 



P = .64 



'13-31 



:: n 



13-31 



:: n 



13-31 



P = .77 



P = .82 

7fT7r.7rj..... 




P = .61 
T7^n . 



■1 1 



2 3-3 -2-10 1 2 3-2 -2 

Continuum Index 



1 1 



Fig. 17. — Histograms comparing the distribution of ri2-Q, ng-i3 and 7113-31 values of objects in 
NGC 1333 (N = 33; shaded) to those in Taurus (N = 58; left), Chamaeleon I (N = 61; middle) 
and L1688 (N = 56; center). P-values from one-dimensional two-sided Kolmogorov-Smirnov test 
are inset in each plot. 
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Fig. 18. — Median spectrum of Class II disks in NGC 1333, alon g with the median s pectra of K5 
M2 objects in Taurus, Chamaeleon I and Ophiuchus (core) from lFurlan et al.l (|2009l ). 
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Fig. 19. — Continuum-subtracted and -normalized 10/im silicate feature profile for the medians 
shown in Figure [18j The profile of pristine dust is scaled to the mean 9.8/im flux of the four regions 
and is shown by the shaded region. 
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Fig. 20.— Dereddened SEDs of TDs and PTDs in NGC 1333 as well as the median of disks in 
NGC 1333, as described in § 15.51 plotted with the dot-dash line. The photosphere for a given 
spectral type scaled to J band is plotted with a dashed line. 



Table 1. Basic Data 



Index 



RA J2000 



Dec J2000 



AORID 



Alternate Namc(s) 



IRAS 03256+3055 

IRAS 2b, [RAC97] VLA 10, SK 7, LAL 157, SVS 19 

SK 14 

ASR32 

IRAS 4b (Southern ouflow cavity) 

IRAS 03262+3123 

LAL 38, [HFR2007] 74 

IRAS 03254+3050 

[HFR2007] 69 

LAL 68 

SVS 9; IRAS 5; ASR 126; LAL 79*, HJ 27, HBC 340, CoKu NGC 1333 3, IRAS 03256+3107, [GMW2007] 9 

ASR67 

ASR 41; LAL 111 

ASR 64; MBO 148 

LAL 166; MBO 38, [GMW2007] 38 

SVS 12; ASR 114; LAL 181; MBO 19, [HFR2007] 45; SK 24; IRAS 6 

SVS 13; ASR 1; LAL 196; HBC 346; SK 13; [HFR2007] 43 

LAL 221; MBO 33 

LAL 223; MBO 51 

LAL 222; MBO 36, [GMW2007] 63 



2 


03 28 45.30 


+31 05 42.1 


22024192 


3 


03 28 57.37 


+31 14 16.2 


22037248 


5 


03 29 04.06 


+31 14 46.7 


22019584 


7 


03 29 11.25 


+31 18 31.7 


22024704 


8 


03 29 12.04 


+31 13 02.0 


22024448 


11 


03 29 23.48 


+31 33 29.4 


22021632 


12 


03 28 32.55 


+31 11 04.8 


22023936 


13 


03 28 34.49 


+31 00 51.2 


22034944 


11 


03 28 34.52 


+31 07 05.5 


22035200 


16 


03 28 39.69 


+31 17 32.0 


22035456 


18 


03 28 43.27 


+31 17 33.1 


22035712 


19 


03 28 48.76 


+31 16 08.9 


22022912 


20 


03 28 51.27 


+31 17 39.5 


22023680 


23 


03 28 57.11 


+31 19 11.9 


22020352 


25 


03 28 58.42 


+31 22 17.7 


22036480 


28 


03 29 01.54 


+31 20 20.7 


22036992 


30 


03 29 03.76 


+31 16 04.0 


22029568 


33 


03 29 08.95 


+31 22 56.3 


22036224 


34 


03 29 09.07 


+31 21 29.1 


22036736 


35 


03 29 09.09 


+31 23 05.7 


22026496 


36 


03 29 10.71 


+31 18 21.2 


22035968 


37 


03 29 12.95 


+31 18 14.6 


22026752 


38 


03 29 18.25 


+31 23 19.9 


22019840 


43 


03 28 38.76 


+31 18 06.7 


22024960 


11 


03 28 43.23 


+31 10 42.7 


22018816 


16 


03 28 44.08 


+31 20 52.9 


22018560 


18 


03 28 47.63 


+31 24 06.1 


22025216 


49 


03 28 47.82 


+31 16 55.3 


22031104 


50 


03 28 51.02 


+31 18 18.5 


22031616 


51 


03 28 51.07 


+31 16 32.6 


22019072 


52 


03 28 51.19 


+31 19 54.9 


22032384 


53 


03 28 52.13 


+31 15 47.2 


22020608 


51 


03 28 52.15 


+31 22 45.4 


22028544 


55 


03 28 52.90 


+31 16 26.6 


22017536 



in 



ASR 30; LAL 261, [GMW2007] 71 

Background Galaxy 

LAL 63 

LAL 80 

LAL 83, [GMW2007] 11 

LAL 93, [GMW2007] 15; MBO 34 

SVS 17; ASR 111; LAL 97; HJ 32, [GMW2007] 16 

SVS 10; ASR 122; LAL 106; MBO 6; Pr 6; LkHa 352A; [GMW2007] 19 

ASR 44; LAL 107; HJ 33; [GMW2007] 20 

ASR 125; LAL 110; MBO 14; [GMW2007] 21; HJ 9; Pr 7 

ASR 45; LAL 125; HJ 42; [GMW2007] 22 

LAL 120; MBO 18; HJ 109; [GMW2007] 23 

ASR 46; LAL 128; HJ 34 



Table 1 — Continued 



Index 



RA J2000 



Dec J2000 



AORID 



Alternate Nainc(s) 



57 


03 28 53.93 


+31 18 09.3 


22031360 


58 


03 28 54.07 


+31 16 54.5 


22025984 


59 


03 28 54.61 


+31 16 51.3 


22030848 


61 


03 28 55.07 


+31 16 28.8 


22030336 


63 


03 28 56.09 


+31 19 08.6 


22029312 


64 


03 28 56.31 


+31 22 28.0 


22029056 


65 


03 28 56.64 


+31 18 35.7 


22032128 


67 


03 28 56.95 


+31 16 22.3 


19056128 


68 


03 28 57.17 


+31 15 34.6 


22023168 


69 


03 28 57.70 


+31 19 48.1 


22028800 


73 


03 28 59.32 


+31 15 48.5 


22030080 


74 


03 28 59.56 


+31 21 46.8 


22032896 


75 


03 29 00.15 


+31 21 09.3 


22022656 


78 


03 29 03.13 


+31 22 38.2 


22033408 


82 


03 29 03.86 


+31 21 48.8 


22033152 


84 


03 29 04.67 


+31 16 59.2 


22026240 


85 


03 29 04.73 


+31 11 35.0 


22021888 


88 


03 29 05.76 


+31 16 39.7 


22030592 


89 


03 29 06.32 


+31 13 46.5 


22027008 


92 


03 29 08.95 


+31 26 24.1 


22018048 


94 


03 29 09.47 


+31 27 21.0 


22020864 


99 


03 29 10.82 


+31 16 42.7 


22027264 


101 


03 29 11.64 


+31 20 37.7 


22023424 


102 


03 29 11.77 


+31 26 09.7 


22019328 


106 


03 29 13.13 


+31 22 52.9 


22033920 


110 


03 29 16.81 


+31 23 25.4 


22021376 


111 


03 29 17.66 


+31 22 45.2 


22033664 


113 


03 29 18.65 


+31 20 17.9 


22021120 


111 


03 29 18.72 


+31 23 25.5 


22034176 


115 


03 29 20.05 


+31 24 07.6 


22034432 


116 


03 29 20.42 


+31 18 34.3 


22031872 


118 


03 29 21.55 


+31 21 10.5 


22025728 


119 


03 29 21.87 


+31 15 36.4 


22029824 


120 


03 29 23.15 


+31 20 30.5 


22027776 



ASR 40; LAL 129; [GMW2007] 25 

SVS 18; ASR 42; LAL 131; [GMW2007] 26; HJ 37 

ASR 43; LAL 136; HJ 38; [GMW2007] 27; [OTS2008] 1 

ASR 107; LAL 141; [GMW2007] 28; HJ 35; [OTS2008] 3 

MBO 146 

LAL 147; MBO 37 

ASR 120; LAL 150, MBO 11; HBC 344; SVS 11; [GMW2007] 30; HJ 21; Pr 8 

SVS 15; ASR 118; LAL 154 

ASR 17; LAL 156; [GMW2007] Sl-5 

ASR 113; MBO 35; HJ 20; [GMW2007] 35 

SVS 16; ASR 106; LAL 171; Pr 11; [GMW2007] 39 

LAL 173; MBO 12; [GMW2007] 41; HJ 8; LkHa 353 

LAL 175; MBO 78; [GMW2007] 42 

LAL 189; [GMW2007] 50; HJ 12; MBO 31 

LAL 195; HJ 11; Pr 14; MBO 5; SVS 8; [GMW2007] 52 

ASR 105; LAL 203; [GMW2007] S2-9 

ASR 99; LAL 205; [GMW2007] 56 

SVS 14; ASR 7; LAL 207; [GMW2007] 58; [RAC97] VLA 22 

ASR 53; LAL 208 

LAL 220; MBO 108 

MBO 60 

ASR 23; LAL 235; [GMW2007] S3A-1 

ASR 85; LAL 245; MBO 58 

LAL 246; MBO 118 

LAL 262; MBO 8; [GMW2007] 72; HJ 15 

LAL 279; MBO 79 

SVS 2; LAL 283; HJ 103; Pr 17; LkHa 270; [GMW2007] 79; MBO 2; HBC 12; VSS 4 

ASR 81; MBO 109 

LAL 293; HJ 102; Pr 18; MBO 16; [GMW2007] 81 

LAL 296; MBO 46; [HFR2007] 67; [GMW2007] 84 

SVS 5; IRAS 03262+3108; ASR 112; LAL 300; [GMW2007] 86; MBO 17 

LAL 304; MBO 32; HJ 17; [GMW2007] 88 

ASR 123; LAL 307; [GMW2007] 90; HJ 46; VSS 27; LkHa 271; SVS 20; HBC 13; StHa 20 

ASR 78; LAL 310; HJ 2; LkH« 355; MBO 30; [GMW2007] 92 



CO 



Table 1 — Continued 



Index 



RA J2000 Dec J2000 AOR ID 



Alternate Name(s) 



121 


03 29 23.23 


+31 26 53.1 


22022144 


122 


03 29 24.08 


+31 19 57.8 


19052544 


123 


03 29 25.92 


+31 26 40.1 


22034688 


125 


03 29 29.79 


+31 21 02.8 


22028032 


126 


03 29 30.39 


+31 19 03.4 


22027520 


127 


03 29 32.56 


+31 24 37.0 


22025472 


128 


03 29 32.86 


+31 27 12.7 


22017792 


131 


03 29 37.72 


+31 22 02.6 


22020096 


133 


03 29 54.03 


+31 20 53.1 


22032640 


134 


03 29 02.69 


+31 19 05.8 


22018304 


135 


03 29 09.41 


+31 14 14.1 


22022400 


136 


03 29 26.79 


+31 26 47.7 


22028288 


137 


03 29 29.26 


+31 18 34.8 


19053312 



LAL 308; MBO 54 

ASR 79; LAL 313; SSTc2d J032924.1+3119 

SVS 4; LAL 318; VSS IX-2; MBO 9; [GMW2007] 98 

LAL 333; HJ 53; MBO 27; [GMW2007] 103 

LAL 336; HJ 49; MBO 28; LkHa 356 

LAL 342; HJ 105; MBO 40; LkHa 357 

LAL 344; MBO 53; [GMW2007] 106 

MBO 65 

ASR 62; LAL 187; MBO 84 

ASR 54 

SVS 6; LAL 321; Pr 19; MBO 15; VSS IX-1; [GMW2007] 100 

LAL 331; SSTc2d J032929.3+3118 



Note. - VLA - [R AC97] - iRodriguez et al.1 
- iHerbig fe Jonesl dl983h . HBC - JHerbig fe Ben] 
lll996h. MBO -IWilking et al.l J2004lV ~j 
' Greissl et"ai] J2007l ^l [OTS20081 



Ijennings et all l ll987l) , SVS 
(J2O10h 



5 



[ SK20011 - ISandell , 
lAspin et al.l Jl994h 



GFT20021 - iGetman et al.l |J2002|), Pr - [P97] 
asa et alj J2008h. 



[HFR2007] - Hatch ell et al 
[SVS76] NGC 1333 - IStrom et all il976h . VSS - IVrba et al 



Kneel [|200ll) , HJ__[HJ83l 
[LAL96] lLada et alj 



03 

o 



LAL 



Preibischl jl997tl . [GMW2007] 
ll2007h . IRAS - [JCC87] 



Jl976ft . SSTc2d - [Merin et al.1 
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Table 2. Spectral Type and Class Information for Sample 



Index 


SpT a 


A v 


Method 6 


Obs ni2-2(> 


Obs n 5 _i2 


State c 


Obs n 2 -25 d 


Class 


2 







neg* 


2.39 


0.12 


envelope 


1.43 


I 


3 








1.26 


-0.62 


envelope 






5 








1.79 


0.23 


envelope 






7 








5.16 


-2.05 


envelope 






11 








2.91 


0.05 


envelope 






12 








1.09 


0.31 


envelope 


0.62 


I 


13 




5.8 


JH 


-0.09 


-0.32 


disk 


0.97 


I 


14 









-0.70 


-0.08 


envelope 






16 




35.7 


HK* 


-0.02 


0.06 


envelope 


0.92 


I 


18 


K7.0 


8.5 


JH 


0.14 


-0.44 


disk 


-0.07 


FS 


19 




17.8 


HK* 


0.20 


0.05 


envelope 


0.38 


I 


20 




5.5 


J 11 


1.44 


0.68 


envelope 


0.33 


I 


2.3 




1.1 


111 


0.78 


0.05 


envelope 


0.22 


FS 


25 




22.2 


JH* 


-1.30 


-0.15 


envelope 


0.62 


I 


28 




10.6 


.JH 


0.72 


0.64 


envelope 


0.89 


I 


30 




13.6 


HK* 


0.25 


0.15 


envelope 


0.23 


FS 


33 




31.6 


HK 


-0.48 


-0.96 


disk 


1.02 


I 


34 




8.1 


.JH 


-1.50 


1.31 


envelope 


0.59 


I 


35 








3.72 


-0.75 


envelope 






36 




2.8 


JH* 


-1.16 


-1.67 


disk 


1.34 


I 


37 








-0.27 


0.48 


envelope 






13 




28.2 


HK* 


-0.02 


-1.05 


disk 


0.05 


FS 


41 




13.3 


JH 


-1.66 


-1.5 


disk 


-1.3 


11 


16 


M7.5 


2.2 


111 


-1.35 


-0.77 


disk 


-1.28 


11 


18 


M4.0 


9.1 


111 


-1.09 


-1.42 


disk 


-1.01 


11 


49 


M8.0 


3 


.JH 


-0.85 


-0.87 


disk 


-0.76 


11 


50 


M4.5 


6.1 


111 


-0.73 


-0.71 


disk 


-0.81 


11 


51 


M4.0 


0.7 


HK 


-1.76 


-1.04 


disk 


-1.32 


11 


52 


K7.0 


5.1 


.JH 


-0.52 


-0.64 


disk 


-0.78 


11 


53 


M7.0 





neg 


-2.62 


-0.9 


disk 


-1.27 


11 


51 


M2.5 


2.7 


HK 


-1.58 


-1.07 


disk 


-1.26 


11 


55 


M5.0 


0.3 


HK 


-1.30 


-1.12 


disk 


-1.16 


11 


57 


K7.0 


14.3 


JH 


0.22 


-0.71 


disk 


-0.49 


11 


58 


M5.0 


1.1 


HK 


-1.72 


-0.51 


disk 


-0.86 


11 


59 


M1.0 


8.1 


.JH 


-0.31 


-0.83 


disk 


-0.65 


11 


61 


M3.0 


10.6 


.JH 


-0.42 


-0.73 


disk 


-0.35 


11 


63 




3.5 


JH* 


-1.43 


-1.72 


disk 


0.99 


I 


61 


M2.0 


14.4 


HK 


-0.24 


-0.24 


disk 


-0.15 


FS 


65 


Ml. 5 


8.1 


111 


0.18 


-0.6 


disk 


-0.7 


11 


67 


M2.5 


9.7 


HK 


0.99 


-0.09 


disk c 






68 


M6.0 


5.4 


HK 


-0.61 


-0.14 


disk c 


-0.42 


11 


69 


M3.5 


3.7 


HK 


-0.45 


-0.13 


disk c 


-0.6 


II 


73 


K2.0 


26 


111 


0.89 


-0.41 


disk 


0.25 


FS 


71 




4.6 


111 


-0.59 


-0.99 


disk 


-0.74 


11 


75 




10.4 


111 


-0.30 


-0.53 


disk 


-0.39 


11 
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Table 2 — Continued 



Index SpT a Ay Method 6 Obs rai2-20 Obs n 5 _i 2 State c Obs n 2 -25 Class 



78 


M2.5 


5.7 


JH 


-0.10 


-0.61 


disk 


-0.24 


FS 


82 


K6.0 


5.7 


JH 


-0.48 


-0.11 


disk c 


-0.58 


11 


81 


M6.0 


8.7 


JH 


0.11 


-0.33 


disk 


-0.05 


FS 


85 




31.7 


HK* 


0.75 


-1.57 


disk 


0.04 


FS 


88 


A3.0 


22.3 


JH 


-1.64 


-0.93 


disk 


-0.63 


11 


89 




43.2 


JH* 


0.02 


-0.98 


disk 


-0.22 


FS 


92 




7 


JH* 


-0.23 


-0.28 


disk 


0.41 


I 


91 




2.1 


HK 


-1.04 


-1.4 


disk 


-0.84 


11 


99 


MO.O 


7.1 


JH 


-0.91 


-0.61 


disk 


-0.13 


FS 


101 


M4.0 


8.1 


HK 


-0.13 


-1.11 


disk 


-0.44 


11 


102 




0.3 


JH 


1.22 


-0.09 


envelope 


-0.02 


FS 


106 


M3.0 


10.1 


JH 


-0.90 


-1.58 


disk 


-0.85 


11 


110 


M2.5 


4.5 


JH 


1.59 


-0.59 


disk 


-0.14 


FS 


111 


MO.O 


3.1 


JH 


-1.22 


-0.96 


disk 


-0.73 


11 


113 




16.1 


JH 


0.74 


0.89 


envelope 


0.32 


I 


111 


MO.O 


3 


JH* 


-0.60 


0.15 


disk c 


-0.4 


11 


115 




26.4 


HK* 


0.39 


0.02 


envelope 


0.6 


1 


116 


K6.0 


13.7 


JH 


0.64 


-0.52 


disk 


-0.07 


FS 


118 


M4.5 





neg 


-0.21 


-0.96 


disk 


-1.04 


11 


119 


K4.0 


4.3 


JH 


-0.97 


-1.31 


disk 


-1.24 


11 


120 


M4.5 


0.2 


HK 


-0.07 


-0.68 


disk 


-0.78 


11 


121 


M4.5 


1.1 


HK 


-1.76 


-1.01 


disk 


-1.02 


11 


122 


K5.0 


6.1 


JH 


2.90 


0.69 


envelope 


0.37 


1 


12.3 




3.8 


JH 


-1.67 


-0.86 


disk 


-1.04 


11 


125 


M4.5 


1.6 


HK 


-0.60 


-0.43 


disk 


-0.91 


11 


126 




1.2 


JH 


0.59 


-0.67 


disk 


-0.69 


11 


127 


M3.0 


7 


JH 


-0.49 


-1.37 


disk 


-0.78 


11 


128 


M4.5 


0.2 


HK 


-0.62 


-1.4 


disk 


-1.15 


11 


131 


M7.0 





neg 


0.93 


-0.88 


disk 


-0.68 


11 


133 


M4.5 


2.1 


JH 


0.33 


-0.45 


disk 


-0.5 


11 


131 




12.2 


JH 


0.53 


-1.99 


disk 


-0.56 


11 


135 








-2.52 


-3.03 


photosphere 






136 


K7.5 


1.8 


HK 


0.39 


-1.93 


disk 


-1.5 


11 


137 


K7.5 


2.96 


HK 


2.66 


-0.51 


disk 


-0.38 


11 



"■Spectral types for #6 8, 84 and 99 fromlGreissl et all J2007h . #110 fro m IWilking et al.1 J2004I) , #135 from 
iGutermuth et al.l J2008bl ), #122 from lMerin et al.l feoirl ), the rest are from I Winston et al.l fcOOSl) 

b Objects with null value in Column 3 had no J, H, or K a band measurements, objects with Ay = 0.00 had 
negative values of Ay, asterix denote objects Ay values which have large uncertainty because their colors are 
outside of the reddened classical T Tauri locus. 

c The state for these objects was determined upon careful examination of their SED and does not correspond 
with the state associated with its ns_i2 value 

d Objects with upper limit for K arc indicated with ... 



63 



Table 3. Continuum indices 



Index 


"2-6 a 


°Vl2-6 


"6-13 


"""e-is 


"-13-31 


CT »H3-31 


18 


-0.72 


0.04 


-0.31 


0.05 


0.08 


0.04 


41 


-2.23 


0.12 


-1.39 


0.14 


-1.32 


0.44 


46 


-1.91 


0.12 


-0.91 


0.18 


-1.16 


0.17 


48 


-1.66 


0.08 


-1.18 


0.14 


-1.00 


0.23 


49 


-1.16 


0.13 


-0.80 


0.15 


-0.78 


0.24 


50 


-1.63 


0.06 


-0.62 


0.07 


-0.91 


0.04 


51 


-1.91 


0.09 


-0.94 


0.14 


-1.08 


0.62 


52 


-1.76 


0.11 


-0.47 


0.14 


-0.89 


0.18 


53 


-1.07 


0.13 


-0.84 


0.15 


-1.93 


0.64 


51 


-1.93 


0.08 


-1.07 


0.12 


-1.12 


0.34 


55 


-1.88 


0.09 


-1.05 


0.18 


0.63 


0.46 


57 


-1.89 


0.11 


-0.50 


0.13 


-0.49 


0.12 


58 


-1.69 


0.14 


-0.52 


0.15 


-0.07 


0.26 


59 


-1.65 


0.06 


-0.58 


0.07 


-0.39 


0.09 


61 


-0.90 


0.05 


-0.59 


0.06 


-0.69 


0.05 


64 


-1.35 


0.19 


-0.17 


0.21 


-0.34 


0.07 


65 


-2.22 


0.10 


-0.33 


0.12 


-0.04 


0.06 


73 


-1.32 


0.08 


-0.21 


0.08 


-0.20 


0.06 


71 


-1.27 


0.06 


-0.97 


0.07 


-0.63 


0.09 


75 


-1.68 


0.19 


-0.81 


0.50 


0.40 


0.48 


78 


-0.65 


0.09 


-0.61 


0.12 


-0.45 


0.18 


84 


-0.82 


0.07 


-0.24 


0.08 


0.38 


0.20 


88 


-1.55 


0.07 


-0.89 


0.09 


-1.29 


0.07 


94 


-1.19 


0.16 


-1.19 


0.67 


-0.64 


0.74 


99 


-0.72 


0.13 


-0.58 


0.17 


-0.88 


0.34 


101 


-1.67 


0.17 


-1.04 


0.36 


1.08 


0.35 


106 


-1.32 


0.06 


-1.46 


0.08 


-0.43 


0.07 


110 


-1.64 


0.23 


-0.94 


0.59 


2.35 


0.55 


111 


-0.86 


0.06 


-0.90 


0.05 


-0.91 


0.06 


116 


-0.95 


0.07 


-0.39 


0.08 


-0.55 


0.08 


118 


-1.88 


0.10 


-0.83 


0.15 


-0.43 


0.26 


119 


-2.01 


0.13 


-1.25 


0.16 


-0.67 


0.13 


120 


-1.56 


0.16 


-0.62 


0.19 


0.02 


0.14 


121 


-1.60 


0.11 


-0.85 


0.13 


-0.77 


0.29 


123 


-1.07 


0.06 


-0.84 


0.09 


-2.19 


0.15 


125 


-1.80 


0.07 


-0.34 


0.09 


-0.92 


0.17 


126 


-1.76 


0.12 


-0.70 


0.18 


0.17 


0.19 


127 


-1.39 


0.13 


-1.30 


0.15 


-0.35 


0.15 


128 


-1.71 


0.08 


-1.34 


0.16 


-0.18 


0.18 


131 


-1.74 


0.09 


-0.76 


0.14 


0.78 


0.14 


133 


-1.29 


1.01 


-0.71 


1.10 


0.26 


0.08 


134 


-1.60 


0.11 


-1.82 


0.12 


0.38 


0.42 


136 


-2.68 


0.08 


-1.73 


0.15 


-0.34 


0.31 


137 


-2.60 


0.07 


-0.22 


0.12 


1.79 


0.09 



^Objects with upper limit for K are indicated with 



Note. — Continuum indices for disk dominated objects in NGC 1333 
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Table 4. Dust Processing Indicators 



Index 


EW 10 a 


<?EW w a 


EW 20 a 


"BW2 ° 


Flo 


CTF 10 b 


F 20 b 


a F 20 b 


Fll.3/F 9 .8 


a Fll.3/F9.8 


18 


1.92 


0.08 


1.52 


0.14 


31.30 


0.89 


14.04 


0.89 


0.53 


0.04 


44 


0.60 


0.06 






0.21 


0.02 


0.27 


0.02 


0.64 


0.14 


46 


1.14 


0.08 


1.16 


0.27 


0.13 


0.01 


0.03 


0.01 


0.77 


0.09 


48 


1.66 


0.08 


1.11 


0.17 


0.89 


0.03 


0.17 


0.03 


0.48 


0.04 


49 


1.07 


0.07 


1.86 


0.18 


1.25 


0.07 


0.59 


0.07 


0.46 


0.06 


50 


4.66 


0.11 


1.88 


0.15 


24.97 


0.31 


3.77 


0.31 


0.52 


0.03 


51 


1.99 


0.09 


2.25 


0.35 


0.23 


0.01 


0.08 


0.01 


0.99 


0.07 


52 


6.05 


0.13 


3.17 


0.17 


14.97 


0.14 


3.07 


0.14 


0.43 


0.02 


53 


1.83 


0.08 






0.55 


0.02 






0.73 


0.06 


54 


3.20 


0.10 


1.35 


0.22 


2.17 


0.04 


0.25 


0.04 


0.55 


0.03 


55 


1.19 


0.08 






0.10 


0.01 






0.54 


0.11 


57 


6.44 


0.14 


3.35 


0.16 


13.23 


0.11 


3.92 


0.11 


0.34 


0.02 


58 


0.81 


0.07 






0.29 


0.02 






1.00 


0.14 


59 


2.94 


0.09 


1.43 


0.15 


6.83 


0.13 


1.53 


0.13 


0.51 


0.03 


61 


2.14 


0.08 


1.88 


0.15 


11.87 


0.31 


4.46 


0.31 


0.50 


0.04 


64 


2.83 


0.09 


2.85 


0.16 


6.75 


0.13 


3.15 


0.13 


0.31 


0.02 


65 


3.13 


0.12 


1.47 


0.15 


7.99 


0.28 


2.53 


0.28 


0.46 


0.03 


73 


8.35 


0.16 


4.87 


0.18 


195.40 


1.26 


77.74 


1.26 


0.28 


0.01 


74 


2.21 


0.08 


1.71 


0.15 


5.56 


0.14 


1.62 


0.14 


0.47 


0.03 


75 


3.37 


0.17 


2.94 


0.29 


0.54 


0.02 


0.30 


0.02 


0.71 


0.06 


78 


3.78 


0.10 


3.87 


0.19 


5.05 


0.08 


2.15 


0.08 


0.48 


0.03 


84 


1.14 


0.07 






0.98 


0.05 






0.36 


0.04 


88 


2.35 


0.08 


1.08 


0.14 


47.59 


1.07 


5.60 


1.07 


0.33 


0.03 


94 


2.31 


0.31 


1.81 


0.34 


0.30 


0.04 


0.10 


0.04 


0.42 


0.06 


99 


1.14 


0.07 


4.42 


0.24 


0.73 


0.04 


0.91 


0.04 


0.44 


0.05 


101 


1.93 


0.11 






0.41 


0.02 






0.46 


0.05 


106 


0.60 


0.06 


1.18 


0.15 


2.19 


0.20 


1.34 


0.20 


0.39 


0.10 


110 


0.74 


0.20 






0.06 


0.02 










111 


0.63 


0.06 


1.00 


0.14 


10.30 


0.91 


4.71 


0.91 


0.49 


0.09 


116 


4.33 


0.22 


4.93 


0.48 


40.16 


1.87 


24.16 


1.87 


0.40 


0.05 


118 


2.36 


0.09 


2.97 


0.25 


0.52 


0.02 


0.24 


0.02 


0.75 


0.06 


119 


1.56 


0.08 


0.60 


0.16 


2.57 


0.10 


0.32 


0.10 


0.56 


0.05 


120 


3.04 


0.10 


3.25 


0.20 


1.42 


0.03 


0.59 


0.03 


0.48 


0.03 


121 


0.73 


0.07 


0.93 


0.27 


0.16 


0.01 


0.04 


0.01 


0.60 


0.08 


123 


4.63 


0.11 


0.44 


0.15 


21.68 


0.26 


1.05 


0.26 


0.55 


0.03 


125 


1.82 


0.08 


0.97 


0.18 


1.02 


0.03 


0.20 


0.03 


0.85 


0.07 


126 


1.12 


0.08 


4.46 


0.25 


0.60 


0.03 


1.02 


0.03 


0.30 


0.05 


127 


2.64 


0.09 


1.77 


0.19 


1.15 


0.03 


0.28 


0.03 


0.45 


0.03 


128 


1.30 


0.09 






0.13 


0.01 






0.97 


0.15 


131 


2.23 


0.09 


0.77 


0.25 


0.15 


0.00 


0.04 


0.00 


0.47 


0.04 


133 


2.26 


0.09 


0.68 


0.15 


2.28 


0.06 


0.58 


0.06 


0.51 


0.04 


134 


2.94 


0.09 


5.51 


0.34 


0.26 


0.01 


0.24 


0.01 


0.45 


0.03 


136 


1.04 


0.08 


7.28 


0.35 


0.43 


0.03 


0.82 


0.03 


0.60 


0.09 


137 


3.08 


0.10 


2.00 


0.15 


0.67 


0.01 


1.16 


0.01 


0.65 


0.04 



'EWj values are in units of /im 



-66- 

b F^ values are in units of 10~ 12 X erg s _1 cm~ 2 

Note. — Objects with unsatisfactory continuum fit are indicated with 
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Table 5. TDs and PTDs in NGC 1333 



Index 



H2-6 vs. ni3_3l 



EWio vs. ni3_3i SED Analysis with Median Classification 



50 


11 


52 


11 


57 


11 


7.3 


II 


101 


PTD 


110 


PTD 


116 


II 


136 


small 712 


137 


TD 



PTD 


PTD 


PTD 


PTD 


PTD 


PTD 


PTD 


PTD 


PTD 


PTD 


PTD 


PTD 


II 


TD 


TD 


TD 


TD 


TD 


PTD 


well above median 


PTD 


11 


well below median 


TD 


TD 


TD 


TD 



Note. — Analysis which indicated the disk was unidentifiable from a radially continuous disk based 
on the given method is denoted with "II" 



